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13.  abstract  (Mdxif^uiti  200  words) 

A  model  of  mass  fire  onset  and  propagation  is  described.  This  is  Version  2  of  a  model  developed  earlier  for 
DNA.  This  computer  model  takes  the  simulation  from  the  detonation  of  one,  or  a  few,  nuclear  weapons  over  an 
urban  site.  The  detonation(s)  ignite  fires,  generally  paper  and  fabrics,  within  rooms  of  houses  and  offices  that 
are  in  view  of  a  burst  point.  Blast  wavc(;>)  extinguish  some  of  the  initial  ignitions,  cause  secondary  ignitions, 
and  can  turn  a  region  from  Ground  Zero  into  rubble.  As  igni.,  '  structures  bum  they,  in  turn,  ignite  other  struc¬ 
tures,  until  an  eventual  lack  of  fuel  allows  the  fire  to  bum  ou'.  The  rate  at  which  fuel  is  used,  i.e.,  the  rate  at 
which  energy  is  released,  governs  the  indraft  surface  winds,  which  also  affect  fire  propagation,  indraft  condi¬ 
tions  can  be  great  enough  that  a  fire  storm  exists.  The  user  describes  an  urban  area  geometrically  as  triangular 
tracks  with  the  number  of  structures  in  each  occupancy  class  (Residential,  High  Rise,  Industrial,  etc.),  locates 
the  weapon  detonation  points,  and  weapon  yields.  Output  requests  produce  graphical  tract  maps  of  wind,  track 
burning  status,  and  energy  release  at  selected  times. 
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SECTION  1 
INTRODUCTION 


This  section  presents  a  brief  history  of  the  Urban  Fire  Simulation  (UFS)  development, 
together  with  a  discussion  of  the  conclusions  drawn  from  the  work  perfonned,  recommendations 
regarding  future  related  work,  and  a  brief  summary  of  this  report. 

1.1  BACKGROUND. 

During  FY  8 1  Mission  Research  Corporation  (MRC)  undertook  to  develop  for  the  Defense 
Nuclear  Agency  (DNA)  a  concept  for  a  model  to  predict  the  onset  and  development  of  mass  fines  in 
urban  areas  as  a  result  of  a  nuclear  attack.  A  review  of  the  literature  showed  that  a  considerable 
amount  of  well  organized  work  had  been  done  toward  the  development  of  a  model  to  predict 
primary  ignitions  resulting  from  nuclear  weapons  detonations,  thk  a  gnat  ^al  of  highly 
fragmented  work  had  been  done  in  studying  phenomena  that  could  be  applicable  to  various  aspects 
of  fire  behavior  dnd  fire  propagation  in  an  utban  environment,  and  that  on-going  work  was 
addressing  flame  extinguishment  due  to  blast  wave  passage,  the  relationship  between  blast  wave 
structural  damage  and  secondary  ignitions,  and  the  relationship  between  building  structural 
characteristics  and  building  burning  characteristics. 

It  was  concluded  that  a  simulation  could  be  developed  m  describe  urban  fire  effects 
resulting  from  a  .;uclear  detonation.  It  was  also  concluded  that  the  highest  risk  components  in  such 
a  simulation  would  probably  be  those  that  dealt  with  the  burning  characteristics  of  buildings,  fire 
effects  on  the  ambient  winds,  and  the  growth  and  propagation  of  uncontrolled  fires  in  an  urban 
environment, 

In  order  to  show  feasibility  of  the  concept  that  was  developed  and  to  provide  a  working 
software  structure  which  would  show  how  the  concept  could  be  implemeatei  a  demonstration 
version  of  the  Urban  Fire  Simulation  (UFS-l)  structure  was  developed  as  part  of  the  FY  81  effort 
(Retercnce  1).  The  demonstration  simulation  was  designed  to  accept  all  initial  ignitions  as  input 
data.  The  UFS-l  development  concentrated  on  the  simulation  structure.  However  development  of 
models  for  the  UFS-l  addressed  the  overall  geometry,  descriptions  of  the  characteristics  of  burning 
buildings  and  burning  blocks  of  buildings  in  a  fire  suppression-free  utban  environment,  effects  of 
the  fire  on  the  ambient  winds,  and  the  effects  of  local  winds  on  the  growth  and  propagation  of  the 
fire. 


The  fire  behavior  and  propagation  models  developed  for  UFS-l  represented  a  first  attempt 
to  put  together  some  to  the  results  of  the  fragmented  fire-related  research  work  into  an  organized 
set  of  models  describing  the  behavior  of  uncontrolled  fires  in  an  urban  environment.  These  models 
were  developed  primarily  for  demonstration  purposes.  As  such,  they  did  not  represent  a  final 
product  and  were  expected  to  be  the  subject  of  careful  review,  revision  and,  where  possible, 
validation. 

During  FY  82  MRC  continued  the  development  for  DNA  of  the  Urban  Fire  Simulation. 
Specifically  the  FY  82  development  addressed: 

•  A  grid  system  to  partition  the  urban  area  into  manageable  pieces  and  methods  of 

incorporating  subgrid  resolution. 

•  Data  to  describe  the  urban  area  for  fire  ignition  and  fire  spread  purposes. 

•  Control  interfaces,  data  interfaces  and  rudimentary  models  for  simulating  ignitions 


1 


related  to  nuclear  weapon  detonation. 

•  Review  and  revision  of  models  simulating  fire  effects  on  the  ambient  winds  and  local 

wind  effects  on  fire  behavior  and  propagation. 

•  Parametric  executions  of  UFS-1  to  show  simulation  sensitivity  to  variations  in  key 

variables. 

The  Uiban  Fire  Simulation  was  expanded  as  part  of  the  FY  82  effort.  The  expanded 
version  (UFS-2)  incorporates,  in  addition  to  the  functions  included  in  UFS-1,  a  more  flexible 
geometry,  together  with  the  data  and  control  interfaces  required  to  support  a  detailed  description 
of  the  urban  area,  primary  ignitions,  blast  extinguishment,  and  secondary  ignitions.  These 
features  have  been  incorporated  into  the  struemre  in  such  a  way  that  a  scenario  can  be  run  starling 
with  weapon  detonation,  proceeding  through  primary  and  secondary  ignitions,  to  fire  spread, 
growth  and  eventual  burnout. 

1.2  CONCLUSIONS  AND  RECOMMENDATIONS. 

The  triangular  grid  used  in  UFS-2  is  much  more  flexible  than  the  rectangular  grid  that  was 
used  in  UFS-1.  The  location,  orientation  and  size  of  each  triangular  tract  can  be  tailor^  to  fit  the 
urban  area  being  analyzed.  Use  of  the  triangular  grid  is  expected  to  alleviate  the  potential  storage 
problem  encountered  in  UFM-1  as  a  result  of  using  uniform  block-sized  rectanguJar  tracts. 

Data  and  control  interfaces  allow  detailed  descriptions  of  urban  area  ignition  and  burning 
characteristics,  given  sufficiently  abundant  data.  Alternatively,  simplified  descriptions  can  be  used 
if  only  sparse  data  are  available. 

The  amount  of  detail  that  can  be  included  in  the  description  of  ignition  characteristics  is 
limited  by  a  combination  of  storage  capacity  and  acceptable  execution  ume.  The  storage  and 
execution  lime  requirements  for  each  traa  are  independent  of  tract  size,  However,  the  total  number 
of  tracts  is  dependent  upon  average  tract  size.  Thus,  execution  time  is  linearly  dependent  upon 
averse  tract  size  and  urban  area  size.  Execution  time  is  strongly  dependent  and  storage 
requirements  are  weakly  dependent  upon  the  number  of  occupancy  classes  in  a  tract.  The 
sensitivity  or  execution  time  to  traa  complexity  can  be  simply  avoided  by  subdividing  complex 
tracts  to  provide  a  larger  number  of  less  complex  tracts. 

The  most  logical  future  directions  for  the  urban  fire  simulation  project  that  are  consistent 
with  th'.  past  development  and  the  need  to  provide  a  product  that  can  be  useful  for  targeting  aie: 

•  The  development  of  appropriately  accurate  and  detailed  models  and  data  to  allow 

prediction  of  primary  ignitions. 

•  Testing  and  validation,  where  possible,  of  the  models  and  the  overall  structure. 

This  is  the  most  conservative  definition  of  the  potential  extent  of  fire  damage.  However,  it  is  also 
probably  the  most  reliable. 

The  development  of  additional  models  and  data  beyond  primary  ignition  to  include  the 
effects  of  blast  damage,  secondary  ignitions  and  fire  spread  is  also  consistent  with  past  work  and 
the  need  to  provide  a  useful  product;  however,  the  potential  utility  of  these  models  and  data  is  less 
clear,  because  of  their  lower  anticipated  prediaive  reliability. 
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1 .3  SUMMARY. 


In  Section  2  the  Urban  Fire  Simulation  (UFS-1)  concept  is  reviewed  and  the  additions  and 
modifications  that  have  been  included  to  define  UFS-2  arc  briefly  described.  These  additions  and 
modifications  affect  the  grid  system,  data  management  system,  subgrid  calculations,  the  events 
used,  and  the  modules  required  to  execute  events. 

In  Section  3  the  triangular  grid  system  is  discussed,  additions  to  the  data  management 
system  are  described,  and  alterations  to  the  Event  Transfer  List  are  identified. 

In  Section  4  the  nuclear  detonation  (NUCDET)  event  is  described,  together  with  the  models 
and  data  involved  in  NUCDET  event  execution. 

In  Section  5  the  blast  effects  (BLSEFF)  event  is  described,  together  with  the  models  and 
data  involved  in  BLSEFF  event  execution. 

In  Section  6  example  model  results  are  presented  and  discussed. 
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SECTION  2 

URBAN  FIRE  SIMULATION  STRUCTURE  AND  CONCEPT 


Ji  this  section  the  UFS-1  concept  is  briefly  reviewed.  The  structural  and  conceptual 
additions  and  modifications  to  UFS-1  are  discussed. 

2.1  BACKGROUND. 

A  detailed  description  of  the  UFS-1  concept  and  structure  is  given  in  Reference  1.  We  will 
present  here  a  brief  summary  to  provide  a  background  for  discussing  the  work  leading  to  UFS-2. 

The  UFS-1  structure  accepted  an  initial  set  of  ignitions  as  input  data  ana  addressed  the 
problem  of  fire  behavior  and  growth  from  that  point.  Since  the  set  of  ignitions  was  a  given, 
ignition  characteristics  were  not  required.  Only  minimal  combustion  characteristics  were  included, 
since  all  executions  considered  a  uniform  fuel  distribution  to  prevent  masking  the  effects  of  other 
variables. 

Each  block-sized  ceil  was  ignited  and  allowed  to  bum  until  burnout  with  a  variable  bum 
rate  that  depended  upon  the  number  of  ignitions  and  the  local  wind  velocities.  Cells  adjacent  to 
burning  ceUs  could  be  ignited  by  convection  or  conduction.  A  fire  was  defined  to  consist  of  a  set 
of  contiguous  burning  cells.  Each  fire  was  subdivided  into  appropriately  symmetrical  subfires. 
Brands  were  considered  to  be  lofted  from  burning  cells  on  the  periphery  of  each  fire. 

The  resoludon  used  in  UFS-1  was  considered  to  be  too  fine.  It  resulted  in  a  limit  on  urban 
area  size  of  about  10,000  blocks  as  a  result  of  requirements  for  storing  cell  data.  The  use  of 
several  ceU  states  required  maintenance  of  a  cell  state  map.  Addidonally,  die  defmidon  and 
maintenance  of  lists  of  fires,  subfircs  and  fire  perimeters  was  required.  The  large  number  of  cells 
used,  together  with  the  considerable  number  of  operations  necessary  to  keep  the  simulation  moving 
forward  cell-by-cell,  and  to  maintain  records  of  the  state  of  each  cell,  fire,  subfirc,  fire  perimeter, 
etc.  involved  the  use  of  a  great  deal  of  computation  resources  and  storage  capacity. 

2.2  CURRENT  STATUS. 

A  flexible  grid  system  has  been  incorporated  into  the  UFS-2  structure.  This  allows  the 
urban  area  to  be  paitidoned  into  triangular  tracts.  The  orientation  of  the  tracts  is  arbitrary,  however 
restriedons  are  placed  on  their  size  and  shape.  Average  tract  size  is  limited  by  the  urban  area  and 
the  maximum  desired  number  of  tracts,  in  view  of  the  related  run  dme  and  storage  requirements. 
Maximum  and  minimum  uact  sizes  are  generally  limited  by  the  variability  of  tract  characterisdes 
and  the  allowed  size  ranges  for  fires,  since  each  tract  is  treated  as  a  fire. 

Nuclear  detonadon  and  blast  effect  events  have  been  incorporated  into  the  UFS-2  structure. 
These  events  include  the  modules  appropriate  to  simulate  primary  and  secondary  ignidons. 

Data  files,  interfaces  and  inidating  roudnes  have  been  included  in  the  UFS-2  stmeture  to 
specify  urban  area  characterisdes  related  to  igniuon  and  combusdon.  Tracts  are  defined  in  terms  of 
fracdonal  occupancy  classes.  Data  are  specified  in  terms  of  "global  data”,  which  apply  to  the  entire 
urban  area  and  "occupancy  class"  data,  which  apply  to  eajh  occupancy  class.  Data  for  each  tract 
are  calculated  from  the  fracdonal  distribution  of  occupancy  classes  in  each  tract  and  the 
distributions  of  characteristics  that  describe  each  occupancy  class.  This  makes  possible  the 
description  of  a  large  number  of  different  tract  types  while  requiring  the  user  to  only  describe  the 
global  and  occupancy  class  characteristics. 
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2.3  TRIANGULAR  TRACTS. 

A  triangular  grid  structure  was  chosen  to  partition  the  urban  area  into  tracts  because  of  its 
flexibility  and  because  it  produces  a  surface  containing  no  unnecessary  step  discontinuities. 
Geometrically,  tracts  of  arbitral'  size,  orientation  and  aspect  ratio  (ratio  of  maximum  to  minimum 
side  length)  may  be  used.  In  this  case  limitations  are  placed  upon  the  tract  size  and  aspect  ratio. 

Symmetrical  fires  are  assumed  by  the  fire  modeling  equations.  Since  each  ignited  tract  is 
treated  as  a  fire,  this  condition  is  satisfied  by  limiting  the  tract  aspect  ratio  to  be  no  greater  than 
two. 


The  average  tract  size  is  determined  by  the  urban  area  and  the  desired  number  of  tracts.  A 
inaximum  tract  area  is  dictated  by  the  minimum  usable  resolution.  Storage  requirements  and  run 
time  are  essentially  independent  of  tract  size.  The  minimum  tract  size  should  be  chosen  to  be 
consistent  with  the  largest  of  the  tract  sizes  dictated  by  storage  requirements,  execution  time,  and 
the  maximum  usable  resolution. 

The  final  number  of  tracts  obtained  will  generally  be  larger  than  the  planned  number.  This 
results  from  subdividing  tracts  to  simplify,  where  required,  the  number  of  occupancy  classes  per 
uact.  It  is  desirable  to  limit  the  number  of  occupancy  classes  per  tract,  because  (NUCDET, 
BSLEFF)  execution  time  for  a  tract  increases  exponentially  with  the  number  of  occupancy  classes 
per  tract,  while  it  increases  linearly  with  the  number  of  tracts  used. 

Figure  1  shows  an  example  of  a  hypothetical  urban  area  that  has  been  partitioned  into 
triangular  tracts.  In  this  case  the  urban  area  had  been  previously  partitioned  into  square  tracts. 
Although  the  resulting  grid  was  used  to  define  the  rectangular  coordinate  system,  the  presence  of 
the  square  tracts  imp^ed  the  definition  of  triangular  tracts  by  making  the  occupancy  class 
boundaries  more  difficult  to  fi'.  with  triangular  tracts. 

Several  attempts,  based  on  different  rationales,  weic  made  to  partition  the  urban  area 
shown  in  Figure  1.  The  ratio.nale  that  served  as  the  basis  for  the  partitioning  shown  in  Figure  1,  and 
used  in  the  example  cdculaiions  is  described  below. 

Define  the  maxiu.um  number  (Nmx)  and  the  minimum  number  (Nmn)  of  tracts  to  be  used. 
Given  the  urban  area  AU.  the  maximum  and  minimum  average  tract  areas  are  given  by 

Amx  =  AU/Ninn  (1) 

and 

Amn  =  AU/Nnu.  (2) 

If  equilateral  tracts  are  assumed,  the  maximum  and  minimum  tract  sides  are  given  by 


and 


i 

Smx  =  (Amx/2)2 
i 

Smn  =  (Anm/2)2. 


(3) 

(4) 


Lay  out  the  tracts  on  the  urban  area  occupancy  class  map  so  mat  tract  boundaries  follow  as  closely 
as  possible  the  boundaries  between  occupancy  classes,  without  violating  (3).  (4)  or  the  maximum 
aspect  ratio  (ARmx).  An  additional  constraint  imposed  on  the  tracts  is  that  each  tract  is  required  to 
have  exactly  three  vertices.  Thus,  it  is  not  possible  for  a  vertex  of  a  tract  to  fall  on  the  side  of  an 
adjacent  tract. 
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Figure  1 .  Layout  and  composition  of  tracts  in  hypothetical  urban  area, 


A  river  flows  from  north  to  south  down  the  west  side  of  the  urban  area  shown  in  Figure  1. 
In  panilioning  the  river  into  triangular  tracts,  the  maximum  aspect  ratio  was  clearly  violated, 
although  (3)  and  (4)  were  observed.  The  aspect  ratio  could  be  violated,  because  these  tracts  cannot 
bum  and.  thus,  will  never  be  represented  as  fires.  It  was  necessary  to  observe  (3)  and  (4),  because 
the  sides  or  the  tracts  defining  the  river  arc  also  sides  of  the  adjacent  tracts,  which  may  contain  fuel 
and  thus  be  represented  as  fires. 

After  the  first  partitioning  the  central  business/  industrial  area  was  found  to  be  contained  in 
a  few  tracts,  each  of  which  exhibited  several  occupancy  classes.  This  area  was  repartitioned  to 
arrive  at  the  configuration  shown  in  Figure  1. 

In  repartitioning  the  central  business/indusuial  area  to  reduce  the  number  of  occupancy 
classes  per  tract,  constraints  (3)  and  (4)  were  related  while  ARmx  was  observed.  Figure  1  still 
shows  tracts  containing  as  many  as  four  occupancy  classes.  These  were  left  in  order  to  be  able  to 
test  the  effect  of  the  number  of  occupancy  classes  on  tract  execution  time. 

Although  triangular  tracts  are  used,  the  computational  coordinate  system  is  still  rectangular. 
Thus,  provisions  must  be  made  to  transform  back  and  forth  between  these  two  coordinate  systems. 
To  accomplish  this,  the  uacts  arc  numbered  as  shown  in  Figure  2  and  the  tract  vertices  are 
numbered  as  shown  in  Figure  3.  The  data  associated  with  tracts  is  contained  in  data  table  TRaCT 
Data  (TRCTD),  and  the  data  associated  with  tract  vertices  is  contained  in  data  table  VERTex  Data 
(VERTD).  The  format  for  data  table  TRCTD  is: 

ITRACT.  IVl,  IV2.  IV3.  NNOC.  lOCT.  FOCT.  TRACT  CENTER.  TRACT  AREA 


where 


rrRACr  is  the  tract  identification  number 

rvi.  IV2.  IV3  arc  the  identification  numbers  of  subject  uact  vertices 

NNOC  is  the  number  of  occupancy  classes  in  subject  tract 

lOCT  is  the  occupancy  class  identification  number 

FOCT  is  the  fraction  of  subject  tract  occupied  by  this  occupancy  class 

TRACT  CENTER  is  the  (X.Y.Z)  coordinates  of  tract  center 

TRACT  AREA  is  the  area  of  tract  in  meters^. 

The  tract  data  table  is  entered  with  the  aact  ID  number  to  obtain  the  above  tract  characteristics.  The 
fomiat  for  data  table  VERTD  is: 

IVERT.  VX.  VY,  VZ.  NAT.  lAT 


where 


rVERT  is  the  vertex  identification  number 

VX.  VY.  VZ  are  the  (X.Y.Z)  coordinates  of  subject  vertex 

NAT  is  the  number  of  tracts  to  which  the  subject  vertex  belongs 

lAT  is  the  set  of  identification  numbers  of  the  tracts  to  which  the  subject  vertex  belongs. 

The  tract  vertex  table  is  entered  with  a  vertex  ID  number  (obtained  from  TRCTD)  lo  obtain  the 
above  vertex  characteristics. 
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Figure  2.  Traci  identification  numbers. 
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2 .4  CONTROL  INTERFACES. 

The  major  modification  that  has  been  made  to  the  control  interfaces  of  UFS-1  is  the  addition 
of  two  new  Events,  the  NUClear  DETonation  (NUCDET)  Event  and  the  BLaSt  EFFects 
(BLSEFF)  Event.  A  number  of  minor  alterations  made  to  the  control  interfaces  will  be  included  in 
the  detailed  descriptions  of  these  events. 

Table  1  shows  the  UFS-2  Event  Transfer  Table,  which  defines  the  intermodule  control 
interfaces.  Hie  NUCDET  and  BLSEFF  Events  have  been  added  to  the  table,  together  with  the 
modules,  models,  data  interfaces,  and  intraraodule  control  interfaces  required  to  execute  the  events. 

2.5  DATA  INTERFACES. 

Three  classes  of  data  are  used  by  UFS-2.  These  arc;  (1)  data  which  arc  user  defined  and 
entered  prior  to  any  execution,  (2)  data  which  arc  calculated  (i.c.  prcproccssed)  from  user  defined 
data  and  entered  prior  to  execution,  and  (3)  data  which  are  calculated  during  execution.  The  data 
sets  belonging  to  Class  1  and  the  approach  used  in  generating  the  data  sets  belonging  to  Class  2 
will  be  described  here.  Descriptions  of  the  data  sets  belonging  to  Class  3  will  be  included  in  the 
descriptions  of  the  modules  and  models  in  which  they  arc  used. 

2.5.1  Class  1  Data  Sets. 

Tract  and  vertex  data,  shown  in  Figure  4,  are  user  specified  from  an  occupancy  class  map 
of  the  urban  area  to  be  aitalyzed.  Currently  these  data  are  manually  specified  and  entered; 
however,  provisions  can  be  made  to  allow  digitizer  entry  of  much  of  this  data  directly  from  the 
occupancy  f'lass  map.  The  above  data  for  each  traa  are  used  by  a  preprocessor  to  calculate  the  area 
of  each  tract  and  the  coordinates  of  the  center  of  each  tract. 

The  dau  which  describe  the  details  of  the  urban  area  for  purposes  of  ignition,  combustion 
and  propagation  of  nuclear  weapon  induced  fires  are  contained  in  a  collection  of  “global”  data  sets. 
The  global  data  sets,  which  are  shown  in  Figure  5,  identify  characteristics  which  apply  over  the 
entire  urban  area.  The  contents  of  the  global  data  sets  arc  currently  u.ser  specified  and  manually 
entered.  However,  the  global  data  are  of  limited  extent  and  they  need  only  be  defined  initially,  and 
updated  when  the  urban  area  characteristics  are  changed. 

The  data  which  describe  the  details  of  each  occupancy  class  in  the  urban  area  for  purposes 
of  ignition,  combustion  and  propagation  of  nuclear  weapon  induced  fires  are  containea  in  a 
collection  of  “occupancy  class"  data  sets.  The  occupancy  class  data  sets,  which  arc  shown  in 
Figure  6,  define  the  distributions  of  characteristics  which  apply  to  each  occupancy  class.  These 
data  sets  are  also  global  in  the  sense  that  the  defined  disuibutions  apply  to  each  occupancy  cl^  no 
matter  where  it  is  found  in  the  urban  area.  The  contents  of  the  occupancy  class  data  sets  are 
currently  user  specified  and  manually  entered.  However,  the  occupancy  class  data  arc  of  limited 
extent  and  need  be  defined  only  for  the  number  of  occupancy  classes  to  be  used  in  the  urban  area  to 
be  analyzed. 

Additional  user  defined  data  include  runtime  data,  grid  data,  and  environmental  data,  which 
arc  shown  in  Figure  7,  together  with  weapon  data,  and  initialization  data  for  fireball  and  irradiance 
calculations,  which  are  shown  in  Figure  8. 

An  effort  has  been  made  to  limit  the  amount  and  detail  of  data  that  the  user  is  required  to 
define,  while  including  enough  descriptive  detail  to  allow  adequate  description  of  the  ignition  and 
burning  characteristics  of  the  urban  area. 
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Table  1 .  Event  transfer  table  for  urban  fixe  simulation. 


READIN 

WRITOT 

SETUP(M) 

UPDATE(M) 

FINISH(M) 

NUCDET 

INPUTS 

OUTPUT 

MODULEOD 

MODULE  (M) 

MODULE(M) 

TPULSE 

INBLST 

WNDMOD 

IGNITE 

BSDEFF 

BSLPRP 

QUENCH 

FULMOD 

SECIGN 

BURN(K) 

ADJCNT(N) 

FUPDT(L) 

STOWV 

STOCS 

STOQ 

BCUPDT(K) 

ACUPDT(N) 

UPDTF(L) 

WVMAP 

eSMAP 

QMAP 

BRNDG(K) 

Define  the  Events  and  Modules  which  appear  in  the  Event-Module  Transfer  Table. 


TERM 


DEFINITION 


READIN 

INPUTS 

WRITOT 

OUTPUT 

SETUP(M) 

MODULE(M) 

UPDATE(M) 

FINISH(M) 

NUCDET 

TPULSE 

INBLST 


WNDMOD 

IGNITE 

BLSEFF 

BLSPRP 

QUENCH 

FULMOD 

SECIGN 

BURN(K) 

BCUPDT(K) 

BRNDG(K) 


Event  to  read  input  data 

Module  to  specify  data  to  be  read  in  and  its  formats 
EveiU  to  write  output  data 
Moc ale  to  specify  data  to  be  output  and  its  formats 
Event  to  setup  module  (M) 

ID  of  module  (M) 

Event  to  bring  the  pertinent  data  elements  of  modulc(M)  up  to  date 
Event  to  teiminate  module(M) 

Event  to  perform  calculations  associated  with  nuclear  weapon  detonation 
Module  to  calculate  free-field  thermal  irradiance  (j/m^)  as  a  function  of  yield, 
burst  height,  and  visibility  at  a  specified  number  of  ranges  from  ground  zero 
Module  to  calculate  free-field  blast  effects  of  peak  overpressure,  positive  phase 
duration,  flow  velocity  and  time  of  arrival  at  a  specified  number  of  ranges  from 
ground  zero  as  a  function  of  yield  and  burst  height 
Module  to  modify  wind  field  to  account  for  column  formation  at  the  burst 
location  as  a  function  of  yield,  burst  height  and  time  after  detonation. 

Module  to  calculate  ignitions  from  this  NUCDET  in  all  affected  tracts 
Event  to  calculate  bto  effects  from  this  NUCDET  in  all  affected  tracts 
Module  to  calculate  free  field  blast  parameters  in  a  specific  tract  and  to  modify 
them  to  account  for  shadowing,  orientation,  and  reflection 
Module  to  calculate  ignitions  extinguished  by  blast  effects  in  a  specific  tract 
Module  to  calculate  blast  induced  foel  modi^ations  in  a  specific  tract 
Module  to  calculate  blast  related  secondary  ignitions  in  a  specific  tract 
Event  to  bring  up  to  date  the  burning  characteristics  of  burning  track  (K) 
Module  to  update  the  data  elements  of  burning  tract  (K) 

Module  to  calculate  brand  generation  and  propagation  for  burning  tract  (K) 
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Table  1.  Event  transfer  table  for  urban  fire  simulation  (continued). 


ADJCNT(N) 

ACUPDT(N) 

FUPDT(L) 

UPDTF(L) 

STOWV 

WVMAP 

STOWCS 

CSMAP 

STOQ 

QMAP 


Event  to  update  the  characteristics  of  tract  (N).  which  is  adjacent  to  a  burning 
tract 

Module  to  update  the  data  elements  of  tract  (N) 

Event  to  update  the  characteristics  of  fire  (L) 

Module  to  update  the  data  elements  of  fire  (L) 

Event  to  store  data  for  a  wind  velocity  map 

Module  to  specify  the  data  formats  for  wind  velocity  maps 

Event  to  store  data  for  a  tract  state  map 

Module  to  specify  the  data  formats  for  tract  state  maps 

Event  to  stom  data  for  a  heat  release  rate  map 

Module  to  specify  the  data  formats  for  heat  release  rate  maps 
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For  each  tract 


Tract  ID  number 

Three  vertex  ID  numbers 

Number  of  occupancy  classes  in  u^t 

ID  number  of  each  occupancy  class  in  tract 

Fraction  of  tract  occupied  by  each  occupancy  class 


For  each  vertex 

Vertex  ID  number 

X,Y.Z  coordinates  of  vertex 

Number  of  tracts  to  which  each  vertex  belongs 

ID  number  of  each  tract  belonging  to  each  vertex 


Figure  4.  User  defined  tract  and  vertex  data. 
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Building  Heights 

Height  of  a  story  (m) 

Number  of  building  heights  used 
Heights  of  buildings  used  (m) 

Building  Classes 

I 

Number  of  building  classes  used 
10  numbers  of  building  classes  used 

Building  Separation  Distances 

Number  of  building  separation  distances  used 
Building  separation  distances  used  (m) 

Foliage-Building  Separation  Distances 

Number  of  foliage-building  separation  distances  used 
Foliage-building  separation  distances  used  (m) 

Foliage  Heights 

Number  of  foliage  heights  used 
Foliage  heights  used  (m) 

Transmissivity  of  each  foliage  height  used 

Occupancy  Classes 

Number  of  occupancy  classes  used 

Building  density  of  occupancy  class  (No.  Bldgsim^) 

Building  fuel  loading  of  occupancy  class  per  floor  (kg/m^/Floor) 
Building  fuel  loading  of  occupancy  class  (kg/m^) 

Room  Furnishings 

Number  of  room  furnishing  types  (critical  ignition  energies)  used 
Critical  ignition  energy  of  each  type  of  room  furnishing  (j^m^) 

Room  Sizes 

Number  of  room  sizes  allowed 
Room  cell  linear  dimension  (m) 

Depth  of  each  '■oom  size  (m) 

Width  of  each  room  size  (m) 


Figure  5.  Global  data  -  user  defined  characteristics  of  the  overall  urban  area. 
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Street  Azimuth  Angles 


Number  of  street  azimuth  angles  used 
Azimuth  angles  of  streets  (rad) 

Street  Widdis 

Number  of  street  widths  used 
Widths  of  streets  (m) 

Window  Sizes 

Window  transmissivity 
Number  of  window  sizes  used 
Height  of  each  window  size  (m) 

Width  of  each  window  size  (m) 

Window  Coverings 

Number  of  window  covering  types  (critical  ignition  energies)  used 
Critical  ignition  energy  of  eich  window  covering  type  (j/m^) 

Figure  5.  Global  data  -  user  defined  characteristics  of  the  overall  urban  area  (continued). 
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Building  Heights 


Number  of  building  heights  in  each  occupancy  class 
ID  number  of  each  building  height  in  each  occupancy  class 
Fraction  of  buildings  in  each  height  class  per  occupancy  class 

Building  Classes 

Number  of  building  classes  in  each  occupancy  class 
ID  number  of  each  building  class  in  each  occupancy  class 
Fraction  of  buildings  in  each  building  class  per  occupancy  class 

Building  Separations 

Number  of  building  separation  distances  in  each  occupancy  class 
ID  number  of  each  building  separation  distance  in  each  occupancy  class 
Fraction  of  buildings  in  each  separation  distance  class  per  occupancy  class 

Foliage-Building  Separation  Distances 

Number  of  foliage-building  separation  distances  in  each  occupancy  class 
ID  number  of  each  foliage-building  separation  distance  in  each  occupancy  class 
Fraction  of  foliage  in  each  separation  distance  class  per  occupancy  class 

Foliage  Heights 

Nurcb-or  of  foliage  heights  in  each  occupancy  class 
ED  nunaber  of  each  foliage  height  in  each  occupancy  class 
Fraction  of  foliage  in  each  height  class  per  occupancy  class 

Room  Furnishings 

Number  of  room  furnishing  t:^es  (critical  ignition  energies)  in  each  occupancy  class 
ID  number  of  each  room  furnishing  type  in  each  occupancy  class 
Fraction  of  room  furnishings  of  each  type  per  occupancy  class 

Room  Sizes 

Number  of  room  sizes  in  each  occupancy  class 
ID  number  of  each  room  size  in  each  occupancy  class 
Fraction  of  rooms  in  each  size  class  per  occupancy  class 

Street  Widths 

Number  of  street  widths  in  each  occupancy  class 
ID  number  of  each  street  width  in  each  occupancy  class 
Fraction  of  streets  in  each  width  class  per  occupancy  class 


Figure  6.  Occupancy  class  data  -  user  defined  characteristics  of  occupancy  classes  over  the  urban 
area. 
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Window  sizes 

Number  of  window  sizes  in  each  occupancy  class 
ID  number  of  each  window  size  in  each  occupancy  class 
Fraction  of  windows  in  each  size  class  per  occupancy  class 

Window  Coverings 

Number  of  window  covering,  tj'pes  (critical  ignition  energies)  in  each  occupancy  class 
ID  number  of  each  window  covering  type  in  each  occupancy  class 
Fraction  of  window  coverings  in  each  type  class  per  occupancy  class 


Figure  6.  Occupancy  class  data  -  user  defined  characteristics  of  occupancy  classes  over  the  urban 
area  (continued). 
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Runtime  Data 


Maximum  simulation  run  time  (s) 
Simulation  tirr.e  to  start  diagnostic  writing 
Number  of  cells  to  be  diagnosed 
ID  numbers  of  cells  to  be  diagnosed 


Grid  Data 

Square  cell  side  length  (m) 

X,Y  coordinates  of  Southwest  corner  of  cell  system 
Number  of  East-West  cells  (columns) 
j  Number  of  North-South  cells  (rows) 


Environmental  Data 

Time  of  data  (s) 

Ambient  wind  speed  (m/s) 

Ambient  wind  direction  (rad) 

Relative  humidity 
Temperature  (deg.  C) 

Low  cloud  layer  height  (m)  above  sea  level 
Low  cloud  layer  type 

High  cloud  layer  height  (m)  above  sea  level 
Surface  albedo 
I  Visibility  index 


Figure  7.  User  defined  mntime,  grid  and  environmental  data. 
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Nuclear  Weapon  Characteristics 

Weapon  (detonation)  ID  number 

X.Y  coordinates  of  detonation  point  (m) 
Height  of  detonation  above  sea  level  (m) 
Total  weapon  yield  (m) 

Detonation  time  (s) 

Maximum  effects  range  (m) 

Minimum  effects  range  (m) 


Radiation  Circle  Rows 

Maximum  number  of  rows  on  a  side  with  no  step 
Maximum  number  of  rows 

Array  of  X.Y  locations  for  the  centers  of  radiation  circle  elements 


Irradiancc  Calculation  Ranges  and  Values 
ID  number  of  weapon 

Number  of  ranges  at  which  iiradiance  is  to  be  calculated 
Value  of  inadiance  at  each  calculation  range 


Figure  8.  User  defined  weapon  data  and  initialization  data  for  fireball  radiation  end  irradiance 
calculations. 
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2.S.2  Class  2  Data  Sets. 


Class  2  data  are  calculated  from  user  defined  data  and  entered  in  the  simulation  data  prior  to 
execution.  These  precalculated  data  are  not  recalculated  for  each  execution,  and,  in  general,  need 
be  recalculated  only  if  the  user  defined  data  arc  altered. 

The  user  defined  tract  and  vertex  data  are  shown  in  Figure  4.  These  data  are  used  to 
calculate  the  area  of  each  tract  and  the  location  of  the  center  of  each  tract.  These  data  are  loaded 
together  with  the  user  defined  tract  and  vertex  data  into  a  single  data  file. 

The  tract  fractional  occupancy  class  specifications  (Figure  4),  the  global  data  characteristics 
(Figure  5),  and  the  occupancy  class  data  characteristics  (Figure  6)  are  us^  in  a  preprocessor  to 
define  the  distributions  of  the  characteristics  for  each  identified  tract  in  the  urban  area. 

For  each  tract  the  fractional  occupancy  class  data  are  of  the  form; 

Ioc(l,  Noc)  =  the  ID  number  of  each  occupancy  class  in  the  tract, 

Foc(Ioc)  =  the  fraction  of  each  occupancy  class  in  the  tract. 

The  data  for  each  global  variable  (sec  Table  2)  are  of  the  form: 

Igv(l,  Ngv)  =  the  ID  number  of  each  value  of  the  global  variable  that  is  used  throughout 
the  urban  area, 

Vgv(Igv)  =  the  magnitude  of  each  value  of  the  global  variable. 

Th ;  data  for  each  occupancy  class  variable  (see  Table  3)  are  of  the  form: 

Iov(l,  Nov)  =  the  ED  number  of  each  value  of  the  global  variable  that  is  used  in  the 
occupancy  class, 

Fov(Iov)  =  the  Fraction  of  each  value  of  the  global  variable  that  is  used  in  the  occupancy 
class. 

Each  occupancy  class  ED  number  (loc)  must  belong  to  the  corresponding  set  of  global  ID  numbers 
(Igv).  Thus,  the  characteristics  that  make  up  any  occupancy  class  data  set  constitute  a  subset  of, 
and  can  be  obtained  from,  the  corresponding  global  data  set.  The  occupancy  class  ID  number 
(lov)  is  used  to  obtain  the  value  of  the  variable  in  ilie  form  Vgv(lov). 

The  number  of  values  of  a  global  variable  in  each  traa  (Ntv)  is  calculated  from  the  number 
of  values  exhibited  by  the  global  variable  in  all  of  the  occupancy  classes  of  the  tract.  The  ID 
numoers  of  the  global  variable  values  in  each  tract  are  the  H)  numbers  of  the  global  variable  values 
in  ail  of  the  occupancy  classes  of  the  tract.  The  fraction  of  each  global  variable  value  in  each  tract 
is  calculated  frxim 


Ftv(Iov) 


Foc(Ioc)  Fov(Iov) 
loc 


(5) 
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where 

Ftv  is  the  fraction  of  the  global  variable  value  in  the  tract, 

Foe  is  the  fraction  of  the  occupancy  class  in  Jie  tract, 

Fov  is  the  fraction  of  the  glob^  variable  value  in  the  occupancy  class. 

The  preprocessor  uses  the  above  calculations  to  define  tract  data  sets  corresponding  to  the 
global  and  occupancy  class  data  sets  which  are  of  the  form; 

It  =  the  tract  ID  number, 

Itv(lJ<tv)  =  die  ID  number  of  each  global  variable  value  in  the  tract, 

Ftv(Itv)  =  the  fraction  of  each  global  variable  value  in  the  tract. 

The  data  sets  which  are  defined  by  the  preprocessor  for  each  tract  are  shown  in  Figure  9. 

2.S.3  Class  3  Data  Sets. 

Class  3  sets  are  calculated  by  the  various  modules  and  models  during  execution.  The 
format  and  contents  of  the  Class  3  data  sets  will  be  described  as  part  of  the  detailed  discussions  of 
the  Events,  mc^ules  and  models  with  which  they  are  associated. 

In  the  following  paragraphs  the  control  and  data  interfaces  defined  for  the  modules  and 
models  associated  with  the  nuclear  detonation  and  blast  effects  events  will  be  discussed.  These 
discussions  will  be  followed  by  a  presentation  and  discussion  of  example  results  from  UFS-2. 
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Building  Height 

Number  of  building  heights  in  tract 
ID  numbers  of  building  heights  in  tract 
Fraction  of  buildings  of  each  height  in  tract 


Building  Classes 

Number  of  building  (structural)  classes  in  tract 
ID  numbers  cf  building  classes  in  tract 
Fraction  of  buildings  of  each  class  in  tract 


Building  Separation  Distances 

Number  of  building  separation  distances  in  tract 
ID  numbers  of  building  separation  distances  in  tract 
Fraction  of  buildings  separated  by  each  separation  class 


Foliage-Building  Separation  Distances 

Number  of  foliage-building  separation  distances  in  tract 
ID  numbers  of  foliage-building  separation  distances  in  tract 
Fraction  of  foliage  separated  by  each  separation  class 


Foliage  Heights 

Number  of  foliage  heights  in  traa 
ID  numbers  of  foliage  heights  in  tract 
Fraction  of  foliage  of  each  height  in  tract 


Room  Furnishings 

Number  of  room  furnishing  types  (critical  ignition  energies)  in  tract 
ID  numbers  of  room  furnishing  types  in  tract 
Fraction  of  room  furnishings  of  each  type  in  tract 


Room  Sizes 

Number  of  room  sizes  in  traa 
ID  numbers  of  room  sizes  in  tract 
Fraction  of  rooms  of  each  size  in  tract 


Figure  9.  Tract  characteristics  -  calculated  characteristics  of  each  in  the  urban  area. 
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Street  Azimuth  Angles 

Number  of  street  azimuths  in  tract 
ID  numbers  of  street  azimuths  in  tract 
Fraction  of  streets  at  each  azimuth  in  tract 


Street  Widths 

Number  of  street  widths  in  tract 
ID  numbers  of  street  widths  in  tract 
Fraction  of  streets  ot  each  width  in  tract 


Window  Sizes 

Number  of  window  sizes  in  tract 
ID  numbers  of  window  sizes  in  tract 
Fraction  of  windows  of  each  size  in  tract 


Window  Coverings 

Number  of  window  covering  types  (critical  ignition  energies)  in  tract 
ID  numbers  of  window  covering  types  in  tract 
Fraction  of  window  coverings  of  each  type  in  tract 


Figure  9.  Tract  characteristics  -  calculated  characteristics  of  each  in  the  urban  area  (continued). 
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SECTION  3 

THE  NUCLEAR  DETONATION  (NUCDET)  EVENT 


The  NUCDET  Event  provides  the  Initialization  for  all  nuclear  weapon  caused  ignitions,  and 
these,  in  turn,  provide  the  initiaLizaiion  for  the  BURN  (fire  behavior  and  spread)  event.  The 
NUCDET  event  includes  modules  for  calculating,  throughout  the  urban  area  of  interest,  the  thermal 
irradiance  (TPUTLSE),  blast  characteristics  (INBLST),  wind  velocity  (WNDMOD),  and  primary 
ignitions  (IGNITE)  resulting  from  detonation  of  a  nuclear  weapon. 

The  stmcture  of  the  NUCDET  event,  and  of  UFS-2  in  general  is  designed  to  accommodate 
single  or  multiple-sequential  nuclear  weapon  detonations.  In  order  for  two  detonations  to  be 
considered  sequential  their  time  spacing  must  be  greater  than  the  arrival  time  of  the  blast  wave 
from  the  fust  detonation  at  the  maximum  range  of  interest.  This  calculational  criterion  decouples 
sequential  detonations  by  not  lowing  the  calculation  of  primary  or  secondary  ignitions  for  the  first 
detonation  to  be  affected  by  the  second  detonation.  Of  course,  both  primary  and  secondary 
ignitions  due  to  the  first  detonation  will  be  ahfected  by  the  second  detonation,  and  these  effects  will 
be  calculated. 

3.1  THERMAL  IRRADIANCE. 

Fire  effects  are  usually  of  interest  over  the  extent  of  an  annulus  centered  on  ground  zero. 
The  inner  radius  of  the  annulus  is  defmed  by  blast  effects  of  such  severity  that  fire  effects  are 
generally  of  little  interest.  The  outer  radius  of  the  annulus  is  defined  by  fire  effects  of  less  severity 
than  specified,  and  essentially  insignificant  blast  effects.  Within  the  annulus  both  blast  and  fire 
effects  are  considered.  This  is  necessary,  because  blast  induced  fuel  modifications  affect 
secondary  ignitions. 

A  potentially  large  number  of  tracts  may  be  included  in  the  annulus,  and  it  is  generally  not 
necessary  to  calculate  the  thermal  irradiance  in  each  traa.  This  is  particularly  mic  if  (as  in  the 
present  case)  it  is  assumed  that  thermal  irradiance  is  a  function  only  of  distance  from  the  burst  point 
G.e.  no  terrain  shadowing).  To  avoid  calculating  thermal  irradiance  for  every  tract  in  the  aimulus, 
thermal  irradiance  can  be  calculated  at  a  user  spewed  number  of  equally  spa^  distances  from 
ground  zero.  Irradiance  data  at  other  distances  from  ground  zero  are  obtained  by  interpolating  on 
these  data.  Alternatively,  a  set  of  monotonic  decreasing  values  of  irradiance  can  be  specified  by  the 
user.  The  corresponding  ranges  from  ground  zero  will  be  calculated. 

Thermal  irradiance  data  arc  calculated  by  module  TPULSE,  which  is  shown  in  Figure  10. 
Upon  entry  to  TPULSE.  environmental  data  (ENVD,  Fig^  7)  are  accessed  to  obtain  the  visible 
range.  Reference  to  Figure  7  shows  that  a  number  of  additional  variables  are  also  included.  Data 
interfaces  have  been  included  to  pemiit  use  of  an  irradiance  model  such  as  described  in  Reference 
2,  although  a  much  simpler  demonstration  model  is  used.  Given  the  visible  range  (Rv),  the 
scattering  coefficient  (Cs)  is  calculated  from  the  definition  of  visible  range  (Reference  3)  by 

Cs  =  0.105/Rv.  (6) 

The  weapon  data  (WEPD,  Figure  8)  and  the  initialization  irradiance  calculation  data  are 
accessed  to  obtain  the  weapon,  detonation,  and  calculation  parameters  needed  to  calculate  thermal 
irradiance  at  specified  locations  in  the  urban  area. 
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Figure  10.  Logic  flow  diagram  for  module  TPULSE. 


cAm.  yncE  #2:  rmn,  pmx,  n'Ques.  and  nques  vaults  of 

IHRADIANCE  (Q(I)  1  ;  GALCLTATE  THE  RANGES  CORRESPONDING  TO 
EACH  OF  THE  Q(I)  v-ALDES 


CALC:  FIREBALL  DATA:  FIREBALL  RADIUS  (RFB),  EFFECTIVE  BURST 
HEIGHT  (HEF) 


CALC:  RADIATION  CIRCLE  CHARACTERISTICS  AS  VIB'ffiD  FROM  EACH  OF 
'CHE  NQUES  CALCULATION  RANGES  FROM  DETONATION;  RADIUS 
(RRCIRd)),  HORIZONTAL  DISTANCE  TO  CENTER  (HDC’)),  VERTICAL 
DISTANCE  'TO  CE^T^E^.  (VD(I)) 


STORE:  IN  RADIATION  DATA  (RRAD):  CALCULATED  IRRADIANCE  DATA  SR(I1, 
KR(I),  Q(I),  RRCIR(I),  HD(I),  VD(I) 


CALC:  TIME  OF  SECOND  THERMAL  MAXIMUM  (T2MX) 


STORE:  'TZMX  IN  WEAPON  DATA  (WEPD) 


Figure  10.  Logic  flow  diagram  for  module  TPULSE  (continued). 


Two  calculation  modes  are  allowed.  In  Mode  1,  at  label  (A)  of  Figure  10,  the  inner  and 
outer  annulus  radii  and  the  number  of  calculation  ranges  are  user  specified.  The  specified  number 
or  equally  spaced  horizontal  (Rh)  and  slant  (Rs)  ranges  from  the  detonation  point,  and  the 
corresponding  free  field  irradiances  at  these  ranges  arc  calculated.  In  Mode  2,  at  label  (B)  of 
Figure  10,  the  inner  and  outer  annulus  radii,  the  number  of  calculation  ranges,  and  the  desired 
inadiances  at  each  of  the  calculation  ranges  are  user  specified.  The  horizontal  and  slant  ranges  are 
calculated  at  which  the  specified  irradiances  will  occur. 

Free  field  thermal  iiradiance  at  each  range  is  calculated  from  the  radiant  exposure-distance 
relationship  (Reference  4,  pp.  316-317) 

Q  =  85.6  (fWt)/(R2)  (7) 

where  Q  is  the  thermal  iiradiance  at  distance  R  (cal/cm^), 
f  is  the  thermal  fraction  of  yield  (taken  to  be  0.3), 

W  is  the  weapon  yield  (KT), 
t  is  the  transmittance, 

R  is  the  distance  from  detonation  point  (Kft). 

At  label  (C)  of  Figure  10  the  fireball  radius  and  radiation  circle  characteristics  are  calculated 
following  (Reference  5,  pp  88-91).  The  fireball  radius  is  given  by 

Rf=K  (YO.35)  Exp(0.465Hb)  (8) 

where  Rf  is  the  fireball  radius  (mi), 

K  is  the  0.53  for  surface  burst, 

K  is  the  0.47  for  Hb  below  40  miles, 

K  is  the  0.41  for  Hb  above  40  miles, 

Y  is  the  weapon  yield  (MT). 

The  radius  of  the  radiation  circle  is  given  by 

I 

Rrc  =  Rf  (1.0-(Rf/Rs)2)2  (9) 

where  Rrc  is  the  radius  of  radiation  circle  (m), 

Rf  is  the  fireball  radius  (m), 

Rs  is  the  distance  to  detonation  point  (m). 

The  horizontal  and  vertical  distances  to  the  center  of  the  radiarion  circle  are  given  by 

1 

Dh  =  Rh  (1.0-(Rf/Rs)2)2  (10) 

I 

Dv  =  Rv  (1.0-(Rf/Rs)2)2  (11) 

Dh  is  the  horizontal  distance  to  center  of  radiation  circle  (m), 

Rh  is  the  horizontal  range  to  center  of  fireball  (m), 

Dv  is  the  vertical  distance  to  center  of  radiation  circle  (m), 

Rv  is  the  venical  range  to  center  of  fireball  (m), 

Rf  is  the  fireball  radius  (m) 

Rs  is  the  range  to  fireball  center  (m). 

These  data  are  stored  in  the  radiation  data  array  RRAD,  which  is  shown  in  Figure  1 1 . 


and 

where 
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Thermal  Inadiance  Data  (RRAD) 


Weapon  (detonation)  ID  Number  (for  each  detonation  in  this  scenario) 

Number  of  ranges  at  which  irradiancc  is  to  be  calculated 
Slant  range  to  detonation  (m) 

Horizonti  range  to  detonation  (m) 

Inadiance  at  this  range  (j/m^) 

Radius  of  radiation  circle  at  this  range  (m) 

Horizontal  distance  to  radiation  circle  center  at  this  range  (m) 
Vertical  distance  to  radiation  circle  center  at  this  range  (m) 

Augmented  Weapon  Data  (WEPD) 

Weapon  (detonation)  ID  number  (for  each  detonation  in  this  scenario) 

X.Y  coordinates  of  detonation  point  (m) 

Height  of  burst  above  sea  level  (m) 

Toti  weapon  yield  (MT) 

Detonation  time  (s) 

Maximum  effects  range  (m) 

Minimum  effects  range  (m) 

Time  of  second  thermal  maximum  (s) 

D  number  of  tract  under  burst  point 


Figure  11.  Thermal  inadiance  and  augmented  weapon  data. 
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The  time  of  second  thermal  maximum  is  calculated  from  the  relationships  (Reference  4,  pp 
309-312) 

T2mxru.C417  wO-‘W,ifHblessthan  15  Kft  (12) 

T2mx  =  0.038  WO'44  (p(h)/p)'3-36,  if  Hb  higher  than  15  Kft  (13) 

where  T2mx  is  the  time  of  second  thermal  maximuni  (s), 

W  is  the  weapon  yield  (KT), 

p(h)/p  is  the  atmospheric  density  ratio  as  a  function  of  height. 

The  time  of  second  thermal  maximum  is  stored  in  the  weapon  data  set  WEPD  and  module 
TPULSE  is  exited.  Figure  8  shows  WEPD  as  it  was  initiated  and  contained  only  user  specified 
data.  Figure  1 1  shows  WEPD  after  execution  of  TPULSE,  when  it  has  been  augmented  by 
inclusion  of  the  calculated  variables  T2mx  and  the  ID  number  of  the  tract  containing  ground  zero. 

3.2  INmAL  BLAST  PARAMETERS. 

Blast  wave  parameters  are  of  interest  for  estimating  the  quenching  effect  of  blast  wave 
passage  on  primaiy  igmtions,  the  interaction  of  the  blast  wave  with  structures,  structure  contents, 
and  primary  ignitions  in  pmducing  secondary  ignitions,  and  the  loss  of  structural  integrity  due  to 
blast  damage. 

Research  into  quenching  of  primary  ignitions  by  blai  wave  passage  (Reference  6)  is  just 
bcginnuig  to  produce  results  that  can  be  used  to  develop  a  predictive  model.  A  predictive  model 
has  been  developed  (Reference  7)  for  secondary  ignitions  resulting  from  blast  wave  ir.teraction 
with  struennes  and  their  contents.  On-going  work  is  expeaed  to  refine  both  the  sfoondary  ignition 
model,  and  the  data  it  uses.  Additional  information  is  required  on  the  effects  of  low  level  blast 
parameters  on  structural  integrity,  and  the  effect  of  loss  of  structural  integrity  on  budding  burning 
charaaeristics. 

Calculation  of  blast  wave  parameters  are  performed  by  module  INitial  BLaST,  which  is 
shown  in  Figure  12.  Blast  wave  parameter  calculations  are  based  on  the  data  prese.nted  in 
Reference  4  on  pages  109  to  121 ,  together  with  the  scaling  ntles  presented  in  Reference  4  on  pages 
100  to  105. 

Upon  entry'  to  module  INBLST,  pressure,  distance,  and  time  scaling  data  are  obtained  from 
the  blast  atmospheric  data  (BLASTD).  BI ASTD  provides  altitude  scaling  factors  at  several 
altitudes  between  sea  level  and  150,000  feet.  These  data  were  obtained  from  Reference  4,  p.  104, 
is  Irivariant,  and  has  been  hard  coded  into  the  module.  Data  are  obtained  from  BLAS'l’D  for 
altitudes  above  and  below  the  altitude  of  interest.  Pressure,  distance  tmd  time  scaling  factors  are 
then  calculated  at  the  altitude  of  interest  by  interpolating  on  these  data. 

Data  file  RRAD  (Figure  1 1)  is  accessed  to  obtain  the  irradiance  calculation  ranges  that  were 
used  in  module  TPULSE.  Blast  wave  parameter.^  are  calculated  at  the  same  set  of  ranges  that  were 
used  for  calculating  irradiance.  As  in  the  case  of  irradiance,  blast  parameters  are  modeled  for 
demonstration  as  functions  only  of  distance  from  ground  zero. 
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GET:  FROM  BLAST  DATA  (BLASTD),  DATA  FOR  CALCULATING  PRESSURE, 
DISTANCE  AND  TIME  SCALING 


CALC:  PRESSURE,  DISTANCE  AND  TI>E  SCALING  FACTORS 


GET:  FROM  RRAD,  IRRADIATION  CAICULATION  RANGES 


Figure  12.  Logic  flow  diagram  for  module  INBLST. 


The  blast  wave  parameters  are  calculated  at  each  irradiance  calculation  range  and  stored  in 
the  blast  data  file  BL5T,  which  is  shown  in  Figure  13. 

3 .3  WIND  MODULE  ( \VNDMOD). 

The  description  of  this  module  is  included  in  the  description  of  the  combustion  and 
propagation  models. 

3 .4  PRIMARY  IGNITIONS  MODULE  (IGNITE). 

The  primary  ignitions  module  IGNITE  uses  the  irradiance  calculated  by  .nodule  TPULSE 
and  the  calculated  tract  characteristics  shown  in  Figure  9  to  calculate  probabilities  of  room  and 
structure  ignition,  and  subsequent  flashover.  Module  IGNITE  is  divided  into  two  parts;  a  part  that 
calculates  global  ignition  (GLBIGN)  variables,  and  a  part  that  calculates  tract  ignition  (TRTIGN) 
variables. 

Ignition  related  calculations  that  apply  to  ail  tracts  and  occupancy  classes  are  performed  by 
submodule  GLBIGN  in  order  to  set  up  global  probability  disuibution  of  the  irradiance  of  window 
coverings  and  of  room  furnishings  as  a  function  of  selected  variables.  Ignition  related  calculations 
that  apply  to  each  tract  and  to  the  occupancy  classes  it  contains  are  performed  by  submodule 
TRTIGN  in  order  to  establish  probability  distributi  ve  of  room  ignitions,  structure  ignitions,  and 
structure  flashovers. 

3.4.1  Global  Ignition  Calculations  (GLBIGN). 

The  thermal  irradiance  of  window  coverings  and  room  furnishings  is  calculated  at  each  of 
the  irradiance  (and  blast)  calculation  ranges  relative  to  ground  zero.  Tnesc  results  will  be  used  in 
TRTIGN  to  determine  the  irradiance  of  window  coverings  and  room  furnishings  in  each  tract, 
taldng  into  account  tract  geometry,  shadowing,  etc.  By  using  the  global  calculations,  the  window 
covering  and  room  furnishing  irr^iance  calculations  ate  only  performed  for  five  to  ten  ranges  from 
ground  zero,  instead  of  for  each  of  potendally  several  hundred  tracts. 

A  simplified  calculation  of  window  covering  and  room  furnishing  irradiance  is  used  for 
demonstration.  The  radiation  source  is  assumed  to  be  on,  and  b«  perpendicular  to,  the  Une 
connecting  the  window  center  and  the  radiation  circle  center,  as  shown  in  Figure  15.  Later,  iu 
TRTIGN,  a  correction  is  made  to  these  results  to  account  for  orientation  of  the  building  with 
respect  to  the  radiation  source.  A  more  detailed  calculation  in  GLBIGN  would  include  a  loop  over 
several  azimuth  angles.  These  data  could  then  be  interpolated  in  TRTIGN  to  obtain  irradiance  as  a 
function  of  buildii.g  orientation. Module  GLBIGN  is  shown  in  Figure  14.  Upon  entry  to  GLBIGN 
the  following  data  are  accessed; 

The  radiation  data  calculated  by  module  TPULSE  and  stored  in 

data  file  RRAD  (Figure  1 1 ). 

The  blast  data  calculated  by  module  INBLST  and  stored  in  data 

file  BEST  (Figure  13). 

The  user  defined  weapon  data  stored  in  data  file  WEPD 

(Figure  11). 
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Weapon  (detonation)  ID  number 

Weapon  detonation  time  (s) 

Number  of  ranges  for  blast  calculations 

Time  of  shock  front  arrival  (s) 

Peak  overpressure  (kg/m^) 

Positive  phase  duration  of  overpressure  (s) 

Horizontal  component  of  dynamic  pressure  (kg/m^) 
Positive  phase  duration  of  dynamic  pressure  (s) 

Figure  13.  Blast  wave  parameter  data  (BLST). 


32 


33 


Figure  14.  Logic  flow  diagram  for  submodule  GLBIGN  (continued). 


END:  LOOP  ON  RADIATION  CIRCLE  ROWS 


END:  LOOT  ON  RANGES  TO  GROUND  ZERO 


Figure  1 4.  Logic  flow  diagram  for  submodule  GLBIGN  (continued). 


Figure  15.  Irrzdiance  geometry  for  room  cells  of  two  nested  roonts. 
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Each  room  size  used  in  the  utban  area  description  is  partitioned  into  elements,  and  the 
number  of  rows  and  columns  of  elements  are  identified  in  the  deepest  and  widest  rooms, 
respectively.  For  calculation  of  room  cell  irradiance,  all  rooms  are  assumed  to  be  nested  about  the 
center  of  a  window  which  faces  the  center  of  the  radiation  circle.  The  window  is  oriented 
perpendicular  to  the  line  connecting  the  center  of  the  window  with  the  center  oi  the  radi'''ion  circle, 
as  shown  in  Figure  15.  The  irradiance  of  each  room  cell  is  calculated  for  one  hypothetical  nom 
whose  width  and  depth  are  equal  to  those  of  the  widest  and  deepest  rooms  used  in  the  urban  area. 

A  loop  is  formed  over  the  irradiance  calculation  ranges  that  were  used  in  TPULSE  to 
calculate  the  free  field  irradiance  data  stored  in  RRAD  (Figure  1 1 ). 

The  radiation  circle  is  partitioned  into  rows  and  columns  by  specifying  the  number  of 
elements  to  appear  in  the  first  row.  The  number  cf  rows  and  number  of  elements  per  row  are  then 
c^culated  from  the  number  of  elements  in  the  fu  •  row  and  the  radiation  circle  u  is,  as  shown  in 
Fi^rc  16.  Partitioaing  the  radiation  circle  allows  inclusion  of  the  effects  of  partial  radiation  circle 
visibility  on  the  probability  of  room  ignition.  The  necessity  for  this  can  be  shown  by  considering 
the  fireballs  from  1 MT  and  10  M'.  bombs  detonated  at  heights  that  maximize  their  2  psi 
overpressure  radii.  At  the  2  psi  overpressure  radius,  the  radiation  circle  will  subtend  about  1 8 
degrees  for  the  1  MT  weapon,  and  about  54  degrees  for  the  10  MT  weapon.  These  are  not  extreme 
condidons,  but  represent  a  range  of  cases  that  might  require  from  slight  to  significant  levels  of 
radiadon  circle  parddoning. 

A  loop  is  formed  over  the  number  o*^  rows  in  the  radiadon  circle.  Assuming  uniform 
radiadon  circle  irradiance,  the  irradiance  of  each  row  is  calculated  from  the  total  number  of 
elements  per  row.  The  irradiatice  of  vertical  window  coverings  is  calculated  at  each  irradiance 
calculadon  range,  and  accumulated  for  each  number  of  radiadon  circle  rows  (from  top  to  bottom  of 
the  radiadon  circle).  These  calculadons  assume  no  shadowing,  normal  incidence  in  the  horizontal 
plane,  and  vertical  incidence  as  defined  by  range  from  ground  zero  and  detonadon  height. 

The  resulting  window  covering  irradiance  data  are  stored  in  the  two  dimensional  array, 

QWC(RGZ,  NROWS) 


where 

RGZ  is  the  range  from  ground  zero,  and 

NROW.S  is  the  number  of  radiation  circle  rows  used. 

Loops  are  formed  over  the  number  of  rows  in  the  deepest,  and  the  number  of  columr  s  in 
the  widest  rooms  used  in  the  urban  area. 

Loops  arc  formed  over  the  number  of  window  sizes  used  in  the  urban  area  and  the  number 
of  radiadon  circle  cells  in  the  row  under  evaluadon. 

For  each  irradiance  calculadon  range  and  each  row  of  the  radiadon  circle,  the  irradiance  of 
vertical  room  furnishings  is  calculated  in  each  room  cell  as  a  function  of  each  window  size.  In  the 
case  of  room  furnishings,  the  irradiance  from  each  radiadon  circle  cell  must  be  examined  because 
of  different  degrees  of  horizontal  shadowing  by  the  different  window  sizes.  Verdcal  window 
shadowing  is  accounted  for  by  radiadon  circle  rows.  ITis  possibility  of  shadowing  is  determined 
by  testing  whether  or  not  the  center  of  each  radiadon  circle  cell  is  visible  from  the  center  of  each 
room  cell. 
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The  room  cell  irradiance  is  stored  in  the  five  dimensional  array, 
QRME(RGZ,  NROWS.  ROWRM,  COLRM.  WNDSIZ) 


where 

ROWRM.  COLRiM  ate  the  coordinates  of  the  room  cell,  and 
WNDSIZ  is  the  window  size. 

The  successive  loops  over  radiation  circle  cells  per  row,  window  sizes,  room  columns, 
room  rows,  radiation  circle  rows,  and  ranges  to  ground  zero  are  closed.  Upon  completion  of  ^ 
the  loops,  module  GLBIGN  is  exited. 

3 .4.2  Tract  Ignition  Calculations  (TRTTGN). 

Toe  tract  ignition  calculation  module  makes  use  of  global  occupancy  class,  and  tract 
characteristics,  together  with  calculated  irradiance  data  to  estimate  the  number  of  buildings  that  will 
be  ignited  in  each  tract. 

Upon  entry  to  submodule  TRTIGN,  shown  in  Figure  17,  the  following  data  are  initialized 
in  the  form  of  data  statements: 

A  data  set  which  specifies  the  probability  of  a  room  fire  occurring  given  ignition  of  window 
coverings  (RFWCI), 

A  data  set  which  specifies  the  probabilitv  of  room  furnishings  being  in  each  room  cell 
(FRFRC), 

A  data  set  which  specifies  the  probability  of  a  room  fire  occurring  given  ignition  of  window 
coverings  (RFRFI), 

A  data  set  specifying  the  probability  of  building  flash  over  given  flash  over  of  at  ieast  one 
room  in  the  building  (FBFO). 

These  data  sets  will  be  described  in  the  following  discussions  of  the  parts  of  the  module 
where  they  are  used.  Weapon  and  detonation  data  arc  obtained  from  data  file  W^D  (Figure  1 1), 
and  thermal  irradiance  data  at  specified  ranges  from  ground  zero  are  obtained  from  data  file  RRAD 
(Figure  11). 

A  loop  is  established  over  all  of  the  tracts  in  the  uiban  area,  and  a  file  is  generated 
containing  all  of  the  tracts  whose  centers  lie  within  the  effects  calculation  annulus.  The  effects 
calculation  annulus  is  defined  by  the  user  specified  minimum  and  maximum  effects  ranges  (see 
WEPD,  Figure  11). 

A  loop  is  formed  over  all  of  the  tracts  in  the  effects  calculation  annulus.  At  label  (P)  the 
tract  center  coordinates  and  tract  occupancy  classes  arc  obtained  from  the  tract  and  vertex  data 
(Figure  4).  The  tract  characteristics  (Figure  9)  are  accessed  to  obtain  distributions  of  the  tract 
ignition  descriptors. 

The  vertical  horizontal  and  slant  ranges  are  calculated  from  the  tract  center  to  the  detonation 
point.  The  horizontal  range  from  iraa  center  to  ground  zero  is  used  to  interpolate  the  irradiance 
data  of  RRAD  (Figure  1 1)  and  the  blast  parameters  of  BLST  (Figure  13)  to  the  center  of  the  subject 
tract. 
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Figure  17.  Logic  flow  diagram  for  submodule  TRTIGN  (continued). 
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L£X3P:  OVER  NUMBER  OF  BUILDING  HEIGHTS  IN  THE  TRACT  AND  IDEirTIFY 
THESE  AS  BUILDINGS  ADJACENT  TO  TARGET  BUILDINGS  .ADJACENT 
BUILOIIGS) 


GET:  FROM  GLCBAL  DATA  AND  TRACT  CHARACTERISTICS,  BUILD^G 
SEPARATION  DISTANCE  AND  FRACTION  CF  BUIIDINGS  IN 
TRACT  SEPARATEE  BY  THIS  DISTANCE 

CAIC:  ELEVATION  ANGLE  FROM  MIDFLOOR  CF  TARCET  BUILDIN3  TO 
BOTTCM  CF  RADIATION  CIRCLE 

CMC;  ELEVATION  ANGLE  FRCM  MIDFLOOR  OF  TARGET  BUILDING  TO 
TOP  OF  ADJACENT  BUILDING 

CMC;  Ntn«ER  CF  ROWS  OF  RADIATION  CTRCIE  THAT' ARE  VISIBLE 
ABOVE  THE  HORIZON 


Figure  17.  Logic  flow  diagram  for  submodule  TRTIGN  (coniinued). 
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Figure  17.  Logic  flow  diagram  for  submodule  TRTIGN  (continued). 


43 


CALC.  CALCULATE  AZiMJTH  ANGLE  BETWEEN  SUBJECT  STPEET  AND 
LINE  FROM  TRACT  CENTER  TO  RADIATION  CIRCLE  CENTER 

CAIC;  PilMtJrH  AND  GECMETRY  CORRECTION  FACTORS  FOR  IRPAD- 
lANCE  OF  WINDOW  CC.’ERINGS  AND  IRRADIANCE  OF  ROOM 
FURNISHtlCS 

CAIL:  IRRADIANCE  OF  WllDOW  COVERINGS  Hi  THE  TRACT  GIVEN 
THE  AZIMUTH  AND  GEXlffiTEX  CORRECTION  FACTOR,  THE 
EFEECrrVE  NU^BER  OF  RADIATION  aPCLE  ROWS  VISIBLE, 
AND  WltOOW  TRANSMISSIVITY 


1  LOOP :  OVER  NUMBER  OF  TYPES  OF  WIMXW  COVERINGS  IN  THE  TRACT 


GET:  FRCM  GLOBAL  DATA  AND  TRACT  CHARACTERISTICS,  CRITTCAL 
IGNITION  ENERGY  OF  OTIDOW  COVERING,  AND  FRACTION  OF 
USING  THIS  WIM3CW  COVERING 

CAIC:  FRACnCN  OF  WINDOW  COVERINGS  TK?T  WILL  IGNITE 
GIVEN  WINDOW  OOVERING  IRRADIANCE 


END:  LJOP  OVER  WINDOW  OOVERING  TYPES 


GET:  FROM  GLOBAL  DATA  AN^  TRACT  CHARACTERISTICS,  NIWEER 
CF  ROVE  AND  OOLUMNS  OF  ROOM  ELE^EN^S  IN  ROOM  SIZE, 
AND  FRACTION  OF  ROOMS  THIS  SIZE  IN  TRACT 

CAIC:  FRACTION  OF  ROa>E  THIS  SIZE  THAT  WILL  HAVE  A  ROOM 
FIRE  GIVEN  IGNITION  OF  ROOM  WINDOW  COVERINGo 
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Figure  17.  Logic  flow  diagram  for  submodule  TRTIGN  (continued). 
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Figure  17.  Logic  flow  diagram  for  submodule  TRTIGN  (continued). 


EUD:  UXSP  OVER  ROOM  ELEMEOTS  IN  IHIS  ROOM  SIZE 


EMD:  UXSP  OVER  WINDOW  SIZES  IN  TRACT 


O 


Figure  17.  Logic  flow  diagram  for  submodule  TRTIGN  (continued). 


CAIC:  ACaJMJIAlE  PROBABILITY  THAT  NO  ROOM  WILL  FLASH 
:7Jm  CM  SUBJECT  FLOOR 


END:  LOOP  OVER  WINDOWS  PER  FIXXR 


CAIC;  ACCUMULATE  PRCBABILITY  THAT  NO  ROOM  WIIL  FLASH 
OVER  IN  SUBJECT  BUILDING 


END:  LOOP  OVER  FLOORS  TO  TOP  OF  BUILDING 


CALC:  PROBABILITY  THAT  AT  LEAST  ONE  ROOM  IN  SUBJECT  BUIIDING 
WILL  FLASH  OVER 

GET:  FROM  TRACT  DATA,  DISTRIBUTIONS  OF  TRACT  CHARACTERISTICS 

CALC:  ACCUMULATE  EXPECTED  FRACTION  CF  BUILDINGS  IN  TRACT  IN 
WHICH  AT  LEAST  ONE  ROOM  WILL  FLASH  OVER 


END:  LOOP  OVER  ROOM  SIZES  IN  IRACT 


END:  LOOP  OVER  STREET  AZIMUTH  ANGLES  IN  TRACT 


Figure  17.  l.ogic  flow  diagram  for  submodule  TRTIGN  (continued). 


© 

END:  LOOP  O'.-."'  .TARGET)  BUILDIIJG  HEIOTTS  IK  TRACT 


Figure  17.  Logic  flow  diagram  for  submodule  TRTIGN  (continued). 
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The  radius  of  the  radiation  circle,  as  viewed  from  the  traa  center,  together  with  the  vertical 
and  horizontal  distances  between  the  tract  and  radiation  circle  centers  are  calculated  from  (9) 
through  (11). 

The  azimuth  and  elevation  angles  from  the  tract  center  to  the  center  of  the  radiation  circle, 
and  the  vertical  angle  from  the  center  to  the  top  of  the  radiation  circle  are  calculated  as  shown  in 
Figure  18.  The  radiation  circle  is  assumed  to  be  normal  to  the  line  connecting  its  center  with  the 
u-act  center.  However,  windows,  window  coverings,  room  furnishings,  etc.  are  not  assumed  to 
be  oriented  normal  to  the  radiation  circle-tract  centerline. 

A  loop  is  formed  over  the  heights  of  the  buildings  in  the  tract.  These  buildings  are 
identified  as  the  buildings  for  which  the  probability  of  ignition  is  to  be  calculated  (i.e.  t^et 
buildings).  The  global  building  height  data  (Figure  5)  are  accessed  using  the  tract  building  height 
data  (Figure  9)  to  quantify  the  tract  Building  heights.  The  tract  building  height  data  includes  the  ID 
numbers  of  the  global  tract  building  heights  that  are  found  in  the  subject  tract.  Given  the  building 
height  ID  number,  the  building  height  is  obtained  from  the  globai  building  height  data.  The 
presence  of  all  tract  descriptors  is  identified  in  the  tract  characteristics  data  (Figure  9)  by  the 
appropriate  giobal  ED  number.  Tract  descriptor  quantitative  data  are  obtained  by  entering  the  global 
dab  with  the  ID  number.  This  procedure  is  repeated  in  accessing  each  of  the  tract  characteristics 
(Figure  5)  data  sets. 

The  number  of  floors  in  the  ^et  building  is  calculated  from  the  building  height  and  the 
specified  floor  height.  One  floor  height  has  been  specified.  In  practice,  a  floor  height  could  be 
specified  for  each  occupancy  class. 

A  loop  is  formed  over  the  floors  of  the  subject  target  building  (height).  At  label  (Q)  the 
midfloor  height  and  the  fraction  of  buildings  in  the  uact  having  this  number  of  floors  is  ctdculated. 
In  practice,  the  fraction  of  buildings  in  the  tract  with  this  number  of  floors  would  be  obtaiued  from 
an  empirically  based  data  table,  l^e  iiradiance  of  room  furnishings  is  calculated  in  a  plane  at  the 
midfloor  height.  The  midfloor  height  is  used  here,  although  the  reference  plane  can  be  placed  at 
any  desired  fraction  of  the  floor  height. 

A  loop  is  formed  over  the  number  of  building  heights  in  the  tract,  and  the  subject  buildings 
are  identified  as  the  buildings  adjacent  to  target  buildings  (adjacent  buildings).  The  adjacent 
building  heights  and  fractions  of  adjacent  buildings  in  the  tract  of  each  height  are  obtained  from  the 
tract  charaaeristics  and  the  globai  data. 

The  elevation  angles  are  calculated  from  midfloor  of  the  target  building  to  the  center  and  to 
the  top  of  the  radiation  circle  as  shown  in  Figure  19.  In  Figure  19  it  is  assumed  tha^  uie  target 
building  is  located  at  the  tract  center.  All  of  the  building-related  calculations  assume  that  the  part  of 
the  structure  for  which  calculations  are  being  made  is  also  located  at  the  tract  center.  The  tract 
calculations,  in  general,  assume  that  the  effects  of  distance  from  the  detonation  point  on  blast  and 
thermal  radiation  parameters  can  be  adequately  accounted  for  by  a  calculation  at  each  tract  center. 
That  is,  the  variation  over  the  dimensions  of  each  tract  is  assumed  to  be  small.  Under  this 
assumption,  distance  'variations  over  structures  within  uacts  can  certainly  be  considered 
insignificant. 

It  may  be  of  interest  at  some  time  to  describe  the  distributions  of  the  blast  and  thermal 
irradiance  parameters  over  each  tract.  This  can  be  easily  done,  given  the  values  of  the  parameters 
at  the  tract  centers  and  the  extent  of  each  tract  in  terms  of  its  vertex  locations. 
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Figure  1 8-  Geometry  of  radiation  circle  from  tract  center. 


A  loop  is  formed  over  the  number  of  building  separation  distances  in  the  tract.  The 
building  separation  distances  and  the  fractions  of  buildings  in  the  tract  separated  by  those  distances 
arc  obtained  from  the  traa  characteristics  and  global  data  sets.  The  elevation  angles  are  calculated 
from  midf.oor  of  the  target  building  to  the  top  of  the  adjacent  building  as  shown  in  Figure  20,  and 
to  the  bonom  of  the  radiation  circle  using  the  geomeuy  of  Fi^re  19.  The  number  of  radiation 
circle  rows  that  are  visible  above  the  horizon  is  calculated  using  the  geometry  of  Figure  19,  the 
vertical  extent  of  the  radiation  circle  that  is  visible,  and  the  number  of  rows  defined  for  the  radiation 
circle 


At  label  (R)  a  loop  is  formed  over  the  number  of  foliage  heights  in  the  traa.  The  foliage 
heights,  transmissivities,  and  the  fractions  of  foliage  in  the  tract  of  those  heights  are  obtained  from 
the  tract  characteristics  (Figure  9)  and  global  (Figure  5)  data  sets. 

A  loop  is  formed  over  the  number  of  fobage-building  separation  distances  in  the  tract.  The 
foliage-building  separation  distances  and  the  fractions  of  foliage  separated  from  buildings  by  those 
distances  are  obtained  from  the  tract  characteristics  (Figure  9)  and  global  (Figure  5)  data  sets. 

The  number  of  rows  of  the  radiation  circle  that  are  shadowed  by  foliage  is  calculated  using 
the  geometry  of  Figure  20  where  the  foliage  is  treated  as  a  non-opaque  adjacent  building.  The 
effective  number  of  radiation  circle  rows  visible  (Nre)  considering  shadowing  by  adjacent 
buildings  and  foliage  is  given  by 

Nre  =  Nrv.(Nsb  +  TfNsf)  (14) 

where  Nrv  is  the  number  of  radiation  circle  rows  visible, 

Nsb  is  the  number  of  radiation  circle  rows  shadowed  by  the  adjacent  building, 

Tf  is  the  foliage  transmissivity, 

Nsf  is  the  number  of  radiation  circle  rows  shadowed  by  foliage. 

A  loop  is  formed  over  the  number  of  street  azimuths  in  the  tract.  The  street  azimuth  angles 
and  the  fraction  of  streets  in  the  tract  at  those  azimuth  angles  is  obtained  from  the  traa 
characteristics  and  global  data  sets.  The  angle  is  calculate  between  the  subject  street  azimuth  and  a 
line  from  tract  center  to  the  radiation  circle  center. 

Correction  factors  are  calculat.  J  for  the  irradiance  of  window  coverings  (CFwc)  and  room 
furnishings  (CFrf)  to  account  for  the  effects  of  building  orientation  and  geomeuy.  Three 
representative  building  shapes  are  shown  in  Figure  21;  (1)  narrow  and  deep,  (2)  wide  and  shallow, 
and  (3)  square.  Building  orientation  is  defined  in  terms  of  the  azimuth  angle  (A)  between  the 
incident  r^iaiion  (lo)  and  the  front  face  of  each  building.  The  building  from  face  is  assumed 
normal  to  the  subject  street.  Each  exterior  room  of  each  building  is  assumed  to  have  a  window  in 
each  exterior  v/all.  Thus,  the  coTCciion  factors  will  be  symmetrical  about  normal  incidence  of 
radiation  on  the  front  face  of  the  building. 
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scfwration  distance  between  adjacent  building  or  foliage  and  target  building 
=  lieight  of  adjacent  building  (or  foliage) 

elevation  angle  frocn  tnidfloor  of  Ixiildi.g  at  tract  center  to  top  of  adjacent  buildincj 
(or  foliage) 


Figure  20.  Ocomctry  of  raclialion  circle  shadowing  by  adjacent  buildings  or  foliage. 


A.  Window  Covering  Irradiance 
Building 
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Figure  21.  Geometry  of  room  and  building  irradiance. 


Figure  21 -A  shows  the  effective  irradiance  of  window  coverings  in  each  room  of  each 
building  type.  Only  one  set  of  window  coverings  per  room  arc  considered.  Thus,  in  those  cases 
where  two  vdndows  are  irradiated  the  maximum  irradiation  is  taken.  Figure  21 -B  shows  the 
effective  irradiance  of  room  furnishings  in  each  room  of  each  building  type.  Room  furnishing 
irradiance  is  taken  as  the  linear  combination  of  the  irmdiances  from  each  window  in  those  cases 
where  two  windows  are  irradiated. 

Table  2-A  shows  the  general  form  of  the  total  window  covering  and  room  furnishing 
iiradiances  given  in  Figure  21.  Also,  tliesc  irradiances  are  tabulated  for  tnrec  different  angles  of 
radiation  incidence.  The  total  irradiances  in  Table  2-A  are  reorganized  in  Table  2-B,  togeilier  with 
the  average  window  covering  and  average  room  furnishing  irr^ances  per  room.  These  are 
averaged  over  building  types  for  each  angle  of  incidence,  and  averaged  over  angles  of  incidence. 

From  Table  2-B  the  correction  factors  for  window  covering  and  room  furnishing 
irradiances  are  given  by 

(  =  0.37  COS(A).  for  COS(A)  .GE.  10  SrN(A). 

CFWC  {  =  0.37  SIN(A),  for  SIN(A)  .GE.  10  COS(A),  (15) 

V  =  0.29  (COS(A)  +  SIN(A)),  otherwise, 

and 

CFrf  =  0.34(COS( A)  -h  SlNf  A)).  ( 1 6) 

At  label  (S)  of  Figure  17  the  irradiance  of  window  coverings  at  the  tract  center  for  the 
effective  number  of  radiation  circle  rows  visible  (QWCIy  is  calculated  by  interpolating  on  the  global 
window  covering  irradiance 

QWC(RGZ,  NROWS), 

and  correcting  the  interpolated  value  for  window  transmissivity,  together  with  geometry  and 
building  type  to  obtain 

QWCI  =  Tw  CFwc  QWe  (17) 

where  Tw  is  the  window  uansmissivity, 

A  loop  is  formed  over  the  number  of  window  covering  types  in  die  tract,  which  are  defined 
in  terms  of  their  critical  ignition  energy  levels.  The  iracuon  of  window  covering  types  that  will 
ignite  (FWCI),  given  irradiance  QWC!  is  calculated.  The  fraction  of  rooms  of  the  subject  size  in 
the  tract  that  wiU  have  a  room  fire,  given  ignition  of  window  coverings  is  calculated  from 
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Table  2.  Building  and  room  irradian;  e  for  example  building  types. 


leti  A  »  B  — -  B  51  A/10 

A.  Total  Irradiance  as  a  Functicrr  of  Buildirrg  Type  and  Geonetry 


Gecm. 


C  »S 
S»C 


1 

1 

C»S 

s»c 

C  51  s 

1 

3 

c»s 
s  »c  . 
c  ~  s 

Total  for  Buildina  Tvrie 


Wimow  Coverings  i  Room  Fumisnings 


5*C  +  2**1 
5.2*C 

7*(C+SV2 


2*C  +  5*S 
3*C 
5.2»S 
7*(C+S)/2 


B.  Average  Room  Irradiance 


Window  Coverinas 


.'ctal  Ave.  Irraa. 

irrad.  per  Roan 


5.2*C 

2.4*C 

3.2*C 


2.4*S 

5.2*S 

3.2*S 


3*(C+S) 
3*(C+S) 
2.5* (C+S) 


0.37*CCCA 


0.37’''SINA 


Roan  Fumisrdnas 


5.2*C, 

2.5*C 

3.3*C 


0.29*(CCSA+SINA), 


2.5*S 

5.2*S 

3.3*S 


7*(C+S) 
7* (C+S ) 
6*(C+-S) 


0.3B*CCSA 


0.38*S24A 


0.34*(CCSArSINA) 
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FRFWCI  B  FWC:  RFWCKRNSIZ,  WNDCOV). 


(18) 


The  data  set  RFWCI  defines  the  probability  of  a  room  fire  given  ignition  of  specified 
window  coverings  and  room  size.  This  data  set  was  initiated  at  the  beginning  of  TRTIGN.  In 
practice  RFWCH  would  be  an  empirically  based  data  table  or  model  of  the  probability  of  room 
Ignition  given  ignition  of  window  coverings. 

A  loop  is  formed  over  the  number  of  room  sizes  in  the  tract.  The  room  size,  the  number  of 
rows  and  columns  of  room  cells  in  the  room  size,  and  the  fraction  of  rooms  of  this  size  in  the  tract 
are  obtained  from  the  tract  characteristics  (Figure  9)  and  global  (Figure  S)  data  sets. 

At  label  (T)  a  loop  is  formed  over  the  number  of  window  sizes  in  the  tract.  The  window 
dimensions  a:id  die  fraction  of  windows  having  those  dimensions  is  obtained  from  the  tract 
chaiactcristics  (Figure  9)  and  global  (Figure  5)  data  sets.  A  single  value  of  transmissivity  is 
defined  for  all  window  sizes. 

A  loop  is  formed  over  the  room  cells  in  the  room  size  by  looping  over  the  number  of  rows 
and  the  number  of  colunms  in  the  room  size.  The  irradiancc  on  vertical  room  furnishings  (QRFI) 
is  calculated  for  each  room  by  interpolating  on  the  global  room  furbishing  irradiance 

QRF(RGZ,  NROWS.  ROWRM.  COLRM,  WNDSIZ). 

using  the  defined  values  of  range  to  ground  zero,  effective  number  of  radiation  circle  rows  visible, 
room  row  and  colunui,  and  window  size.  The  interpolated  value  is  corrected  f'^r  window 
transmissivity,  together  with  geometry  and  building  type  to  obtain 

QRH  =  Tw  CFrf  QRF.  (19) 

A  loop  is  formed  over  the  room  furnishings  in  the  tract,  which  arc  defined  in  terms  of  their 
critical  ignition  energy  levels,  The  fracUon  of  room  furnishings  th.it  will  ignite  (FRFI)  given 
inadiance  QRFI  is  calculated.  The  fraction  of  room  furnishings  expected  to  be  found  in  the  subject 
room  cell  is  obtitined  from  daui  set 

FRFRC(ROWRM.  COLRM,. 

which  w.-s  one  of  the  data  sets  initialized  at  the  beginning  of  Submodule  TRTIGN.  In  practice 
FRFRC  would  also  be  a  function  of  room  size,  and  would  represent  an  empirical  average  of  the 
fraction  of  furniture  cue  would  expect  to  find  in  a  given  ceil  of  a  room  of  given  size  in  the  urban 
■irca  being  cvolimtcd. 

The  fraction  of  room  famishing'^  that  will  ignite(FRMFI),  given  irradiance  QRFI  and  room 
furnishings  present  is  accumulated  by 

y  y 

s  ^  ^  FRR  FPFRC  (RMROW,  RMCOL)  (20) 

RMROW  RMCOL 

tur  all  of  the  room  cell-.  In  the  subject  room  size. 
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At  label  (U)  the  loops  on  room  columns  and  rows  (i.e.  room  cells),  and  window  sizes  are 
terminated. 


The  fraction  of  rooms  in  the  tract  of  the  subject  size  that  will  have  a  room  fire  given  ignition 
of  room  furnishings  (FRFRPI)  is  calculated  in  terms  of  the  fraction  of  room  furnishings  ignited 
(FRMFI)  and  the  probability  of  a  room  fire  in  a  room  of  this  size  given  ignition  of  room 
furnishings.  Thus, 

FRFRH  =  FRMH  RFRHCRMSIZ,  FRH).  (21 ) 

The  data  set  RFRFI  was  initiated  at  the  beginning  of  TRTIGN.  In  practice  this  would  be 
an  empirical  data  set  or  model  based  on  experimental  or  experiential  data  defining  the  relationship 
between  room  furt^hing  ignition  and  occurrence  of  a  room  fire. 

A  loop  is  formed  over  the  number  of  building  classes  in  the  tract.  The  building  class  and 
the  fraction  of  buildings  in  the  tract  or  this  class  are  obtained  from  the  tract  characteristics  (Figure 
9)  and  glolwil  (Figure  5)  data  sets. 

The  fraction  of  rooms  that  will  flash  over,  given  a  room  fire  due  either  to  ignition  of 
window  coverings  or  to  ignition  of  room  furnishings  is  calculated  from 


FFORF  =  FBCL  FBFo  (1  -  0.5  ((1  -  FRFWCI)  + 

(1  -  FRFRFI))) 

where  FBCL  is  the  fraction  of  buildings  in  the  tract  of  this  class, 
FBFD  is  the  fraction  of  buildings  in  the  tract  of  this  class 
that  will  flash  over  given  at  least  one  room  flashed  ovei. 


(22) 


A  loop  is  formed  over  the  number  of  floors  from  the  subjea  floor  to  the  top  of  the  building. 
In  this  case  the  distribution  of  windows  in  the  tract  per  floor  (DAWF(FLOORHT))  is  calculated 
under  the  assumption  that  the  number  of  buildings  of  given  height  decays  exponentially  with 
increasing  building  height.  In  practice  this  would  be  an  empirically  derived  distribution  for  the 
given  urban  area. 


The  probability  that  no  room  on  the  subject  floor  will  flash  over  is  calculated  at  label  (V) 

from 


PNRFOF  =  ^  DWPF(FLOORHT)  (1  -  FFORF)NWPF(FLOORHT) 
FLOORKT 


(23) 


where  NWPF  is  the  number  of  windows  on  the  subject  floor.  The  loop  over  floors  to  the  top  of 
the  building  is  terminated  and  the  probability  that  no  room  in  the  subject  building  will  flash  over  is 
calculated  from 


PNRFOB  = 


NTBF 


PNRFOF. 


(24) 
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The  probability  that  at  least  one  room  in  the  building  will  flash  over  is  given  by 

PIRFOB  =  1  -  PNRFOB.  (25) 


The  fraction  of  buildings  in  the  traa  that  will  flash  over  is  accumulated  and  ncrmalizcd  by 
including  the  distributions  of  the  tract  chaiacteiistics  and  summing  t24)  over  all  of  the  defined  tract 
characteristics.  The  complete  expression  is  given  by 


FBFOT  = 


NRMS  NS.\Z  NFSP  NFHT  NBSP  NABH  NTBF  NTBH 


(26) 


(FRMS  FSAZ  FFSEP  FFOLHT  FBSP  FB2frr  FBIHT  F31FHT) 


(1  ■ 


j  J^TBKTI  ( 


TBKTI  (  ^WPF(FWPF  (1  -  PIRFOB))). 
NTBF  NWPF 


The  loops  are  terminated  over  tract  room  sizes,  street  azimu'.n  angles,  foliage-building 
separation  distances,  foliage  heights,  building  separation  distances,  adjacent  building  heights, 
number  of  floors  in  target  buildings,  and  target  building  heights. 

The  number  of  new  buildings  flashed  over  in  the  subject  tract  is  calculated  from 

NNEW  =  FBFOT  NBT  (27) 

where  NBT  =  the  number  of  buildings  in  the  subject  tract.  The  total  number  of  buildings  flashed 
over  in  the  tract  (NTOT)  is  updated  by  accumulating  NNEW  to  NTOT.  The  number  of  buildings 
that  have  flashed  over  during  the  most  recent  171  seconds  (NLAST)  is  updated  to  the  current  time 
(T2MX)  by  accumulating  NNEW  and  deleting  those  that  no  longer  satisfy  the  criterion. 

The  new  values  of  NTOT,  NNEW,  NLAST,  and  TLAST  are  stored  in  the  ignition  data 
(IGND,  Figure  22),  and  submodule  TRTIGN  and  module  IGNITE  are  exited. 


3.5  THE  BLAST  PROPAGATION  MODULE  (BLSPRP). 

The  blast  propagation  module  specifics  the  ulast  effects  event  parameters  for  each  tract 
within  the  annulus  defined  by  RMXEFF  and  RMNEFF.  The  demonstration  module  simply 
interpolates  on  the  global  BLST  data,  assuming  free  field  conditions.  In  practice  ihe  module  could 
take  into  consideration  the  effects  of  terrain,  structures,  and  forests  on  blast  wave  parameters. 

Upon  entry  to  the  BLSPRP  module,  as  shown  in  Figure  23,  weapon  data  are  accessed  to 
obtain  coordinates  of  the  detonation  point.  These  are  ased  to  identify  the  tract  that  contains  ground 
zero. 

The  global  irradiance  and  blast  data  are  accessed,  and  a  loop  is  formed  over  the  tracts  in  the 
annulus  bounded  by  the  user  specified  minimum  and  maximum  effects  ranges.  The  file  containing 
these  tracts  was  created  in  TRTIGN.  and  can  be  accessed  by  BLSPRP.  In  general,  any  data  file 
created  by  a  module  in  an  event  can  be  accessed  by  all  subsequent  modules  in  the  event. 
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Tract  ID  number 

Time  of  first  ignition  in  tract 

Total  number  of  ignitions  in  tract 

Number  of  ignition  times  during  the  last  171  seconds 

Times  of  ignitions  during  the  last  171  seconds 

Number  of  ignitions  at  each  ignition  time  during  the  last  171  seconds 


Figure  22.  Tract  ignition  data  (IGND). 
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Figure  23.  Logic  flow  diagram  for  module  BLSPRP. 
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The  global  thermal  irradlonce  and  blast  parameters  ore  interpolated  to  the  range  of  the 
subject  tract  center  from  ground  zero.  Logarithmic  interpolation  is  used,  assunting  square  law 
propagation  of  thermal  irradiance  and  blast  wave  effects,  together  with  free  field  conditions.  The 
interpolated  values; 

TASF  a  lime  of  shock  front  arrival, 

POP  =  peak  overpressure, 

PPDOP  s  positive  phase  duration  of  peak  overpressure, 

PHDP  s  peak  horizontal  dynamic  pressure, 

PPDHP  =  positive  phase  duration  of  horizontal  pressure, 

QI  =t  thermal  irradiance, 

are  stored  in  the  calculated  tract  data  file  (CATD),  shown  in  Figure  24. 

A  blast  effects  event  for  the  subject  tract  is  placed  on  the  event  list,  shown  in  Figure  24,  to 
be  processed  at  the  interpolated  time  of  shock  front  arrival  at  the  tract  center.  The  loop  over  tracts 
in  the  annulus  is  terminated,  and  module  BLSPRP  is  exited. 

Tills  lerniinates  the  discussion  of  the  NUCDET  event.  The  blast  effects  event  BLSEFF  is 
discussed  in  the  next  section. 
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Calculated  Tract  Data  File  (CATD) 


ID  number  of  tract 

Time  of  shock  front  arrival  at  tract  center  (s) 

Peak  oveipressure  at  tract  center  (g/tifi) 

Peak  horizontal  dynamic  pressure  at  tract  center  (g/'m^) 

Positive  phase  duration,  overpressure  (s) 

Positive  phase  duration,  dynamic  pressure  (s) 

Thermal  inadiance  at  tract  center  (jAn^) 

Event  Data  Hie  (EVNT)  Format  for  BLSPRP  Event 

Simulation  time  at  which  event  is  to  be  processed  (s) 

Event  ID  number 

Weapon  (detonation)  ID  number 

Tract  ID  number 

Figure  24.  Calculated  tract  data  (CATD)  and  event  data  (EVNT). 
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SECTION  4 

THE  BLAST  EFFECTS  (BLSEFF)  EVENT 

The  blast  effects  event  addresses  the  quenching  of  existing  ignitions,  structural  damage, 
and  secondary  ignitions  resulting  from  passage  of  the  blast  wave. 

The  effect  of  blast  wave  passage  on  existing  ignitions  is  estimated  in  module  QUENCH 
and  the  effect  on  structural  integrity  is  estimated  by  module  FULMOD.  Secondary  ignitions 
resulting  from  the  interaction  of  blast  parameters  with  stmctures  and  existing  fires  is  estimated  by 
module  SECIGN. 

The  objective  of  the  blast  effects  module  is  to  specify  the  total  number  of  primary  and 
secondary  ignitions  resulting  from  detonation  of  a  nuclear  weapon  in  an  urban  area.  The  structure 
of  the  blast  effects  event  is  deigned  to  accomraodme  sequential  detonations  as  defined  in  Section 
3.0. 

4.1  CALCULATION  OF  BLAST-FIRE  INTERACTIONS  (QUENCH). 

Blast-fiie  interactions  have  been  the  subject  of  considerable  theoretical  and  empirical 
Investigation  (References  7  and  8).  These  investigations  are  just  beginning  to  yield  results  upon 
which  a  model  can  be  based.  The  procedure  used  in  module  QUENCH  is  based  on  the  work  of 
Martin  reported  in  Reference  8. 

Module  Quench  is  shown  in  Figure  25.  Upon  entry,  the  ignition  data  set  IGND  (Figure 
22)  is  accessed  to  obtain  tract  ignition  data,  and  the  calculated  tract  data  set  CATD  (Figure  24)  is 
accessed  to  obtain  tract  blast  wave  parameters. 

Only  ignitions  that  have  occurred  in  the  past  171  seconds  are  considered  subject  to  being 
quenched  by  blast  wave  passage.  Thus,  if  there  arc  no  ignition  times  in  the  past  171  seconds,  the 
module  is  exited. 


If  ignition  times  arc  found  in  IGND  during  me  past  171  seconds,  the  probability  is 
calculated  that  they  will  be  quenched  by  the  blast  wave.  The  calculation  of  the  probability  of 
quenching  an  ignition  is  based  on  Figure  26. 

The  boundaries  (B1  and  B2)  between  die  "quench"  and  "no  quench"  regions  of  Figure  26 
are  calculated  from  the  time  since  ignition  by 


B1  =  0.0486  (TASF  -  TLAST(I)) 

B2  =  14.02  -  0.0777  (TASF  -  TLAST(I)) 


(28) 


where  TASF  is  the  shock  front  arrival  lime,  and 

TLAST(l)  is  the  1-lh  ignition  time  in  the  past  171  seconds. 


1  GET:  FROM  IGNITTCN  DATA  SErT  IGND,  IGNITION  DATA  _  I 

i  ;  FROM  OMCULATED  TRACT  DATA  SET  CATD,  TEIRCT  BLAST  V3AVE  PARAMETERS  i 


Figure  25.  Logic  flow  diagram  for  module  QUENCH. 
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Based  on  the  peak  overpressure  experienced,  the  probability  of  the  ignition  being 
extinguished  is  calculated  from 


Pex  = 


0.25  if  POP  .LT.  5.4  (psi)  .AND.  POP  .LT.  Bl, 

0.78  if  POP  .LT.  5.4  (psi)  .AND.  POP  .GT.  Bl, 

0.94  if  POP  .GT.  5.4  (psi)  .AND.  POP  .LT.  B2, 

0.15  if  POP  .GT.  5.4  (psi)  .AND.  POP  .GT.  B2. 

The  number  of  ignitions  at  TLAST(I)  that  will  be  quenched  by  the  blast  wave  is  calculated  from 
Nex  =  Pex  NLAST(I).  (31) 


(29) 


(30) 


The  loop  over  late  ignition  times  is  terminated,  the  values  of  total  (NTOT)  and  late 
(NLAST)  ignitions  are  modified  by  Nex.  the  new  values  are  stored  in  the  ignition  data  set  IGND, 
and  module  QUKnCH  is  exited. 

4.2  BLAST  EFFECTS  ON  FUEL  PROPERTIES  OF  STRUCTURES  (FULMOD). 


The  modification  of  smictures  by  the  passage  of  the  blast  wave  is  of  interest  because  it 
contributes  to  secondary  ignitions  and  it  affects  the  way  a  structure  will  bum.  Blast  wave  damage 
results  in  loss  of  structural  integrity,  which  in  turn  results  in  increased  fire  ventilation  which  can 
reduce  the  fraction  of  rooms  that  may  flash  over.  This  causes  a  slower  building  burning  rate  and 
lower  burning  temperatures,  with  reduced  impact  on  the  ambient  air  movement. 

Structural  fire  damage  is  defined  as  a  function  of  peak  overpressure  at  the  tract  center 
following  the  procedure  described  in  Reference  9,  pp.  38-47,  The  results  of  structural  damage 
studies  and  tests  that  are  available  and  on-going  (Inferences  6, 10, 11)  can  be  used  to  improve  this 
model. 


Upon  entry  to  module  FULMOD,  shown  in  Figure  27,  the  tract  building  classes  are 
obtained  from  the  tract  characteristics  (Figure  9),  and  the  blast  wave  parameieis  at  the  tract  center 
are  obtained  fi'om  the  calculated  tract  data  (Figure  24) 

A  loop  is  formed  over  the  number  of  structural  building  classes  in  the  tract.  Ten  building 
classes  are  included,  these  are  defined  in  Reference  9,  pp.  25-32.  Within  the  loop  over  building 
classes,  a  loop  is  fonned  over  the  number  of  struCTure  damage  classes  used.  Currently  three 
damage  classes  are  used,  these  are  defined  in  Reference  9,  pp.  38-46.  Given  the  peak 
overpressure,  the  damage  level  for  each  building  class  is  calculated  following  the  procedure  shown 
in  Reference  9,  p.  47. 

The  structural  damage  level  for  each  budding  class  in  the  tract  is  stored  in  the  calculated 
tract  data  (CATD,  Figure  24),  and  module  FULMOD  is  exited. 

4.3  CALCULATION  OF  SECONDARY  IGNITIONS  (SECIGI.;. 

The  calculation  of  secondary  ignitions  is  based  on  the  procedure  described  in  Reference  9, 
pp.  49-64.  Reference  9  employs  an  alignment  chart  that  provides  an  estimate  of  the  number  of 
secondary  ignitions  that  will  occur  in  a  given  building  structural  class  as  a  function  of  building 
contents  and  blast  induced  damage  level.  The  alignment  chan  has  been  represented  in  the 
simulation  by  a  set  or  exponentii  functions. 
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QTIER 


GET:  mXA  TRACT  CHARACTERISTICS,  BUILDING  CLASSES  IN  THIS  TRACT 
:  FROM  CALCULATE)  TRACT  DATA  FILE  (CATD),  BLAST  WAVE 
PARAMETERS  AT  TRACT  CENTER 


Figure  27.  Logic  flow  diagram  for  module  FULMOD. 


upon  ena7  to  module  SECIGN,  shown  in  Figure  28.  the  tract  characteristics  data  is 
accessed  to  obtain  the  number  of  building  classes  and  building  contents  in  the  tract,  together  with 
their  distributions.  The  calculated  tract  data  (CATD)  is  accessed  to  obtain  the  damage  level 
experienced  by  each  structural  building  class.  The  ignition  data  (IGND)  is  accessed  to  obtain  the 
number  of  total  and  late  ignitions  in  the  tract. 

A  loop  is  formed  over  the  number  of  structural  building  classes  in  the  tract.  The  total 
number  of  secondary  ignitions  (NSEC)  is  calculated  from 

NSEC  =  X  AO  EXP(A1  ITBCL)  EXP(A2  ITBCN)  NBT  FTBCL  (32) 

ITBCL 

where  AO,  A1  and  A2  is  the  coefficients  defined  by  the  damage  level. 

ITBCL  is  the  smictural  building  class  ED  number, 

FTBCN  is  the  building  contents  ID  number, 

NBT  is  the  total  number  of  buildings  in  tract, 

FTBCL  is  the  fraction  of  buildings  in  tract  of  subject  smictural  building  class. 

The  loop  over  building  classes  is  terminated,  the  numbers  of  total  and  late  ignitions  in  the 
tract  are  modified  in  accordance  with  NSEC,  the  new  values  are  stored  in  the  ignition  data,  and 
module  SECIGN  is  exited. 

Module  SECIGN  is  the  final  module  in  the  execution  of  the  BLSEFF  event.  At  this  pouit 
the  number  of  ignitions  stored  in  IGND  represents  the  total  number  of  ignitions  in  each  traa  in  the 
annulus  bounded  by  RMNEFF  and  RMXEFF  that  were  caused  by  detonation  of  the  subject 
nuclear  device. 
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eitter 


GET:  FRCM  TRACT  CHARACTERISTICS,.  TRACT  BUILDING  CLASSES 
AND  BUILDING  CONTENTS 


GET:  FROM  CAICULATED  TRACT  DATA  (CATD) ,  DAMAGE  LEVEL  FOR  EACH 
BUILDING  CLASS 


GET:  FROM  IGNITICN  DATA  (IGND),  NtMBER  OF  IGNITIONS  IN  TRACT 


LOOP;  OVER  NUMBER  OF  BUIIDING  CLASSES  IN  TRACT 


CALC:  NUMBER  OF  SECONDARY  FIRES  IN  BUILDING  CLASS 
GIVEN  BUILDING  CONTENTS  AND  DAMAGE  LEVEL 


Figure  28.  Logic  flow  diagram  for  module  SECIGN. 
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SECTION  5 


RESULTS  FROM  I'lUCLEAR  DET^  NATION 
AND  BLAST  EFFECTS  EVENTS 


The  nuclear  detonailon  and  blast  effects  events  have  been  executed  many  times  during  the 
debugging  prtxess.  The  results  of  some  of  these  executions  are  presented  here  to  show  the  types 
of  results  that  these  events  wUl  provide. 

Figures  29, 30,  and  31  present  results  for  tract  76  from  three  executions  of  the  NUCDET 
and  BLSEFF  events  for  three  values  ot  weapon  yield  and  detonatioti  altitude.  Tract  76  is  100 
percent  residential  occupancy.  7.5  million  souare  meters  in  area,  and  contains  5002  structures. 
Figure  29  shows  thermal  irradlancc  vetms  distance  from  ground  aero  as  a  function  of  weapon 
yield  and  dettntatlcai  altitude.  It  can  be  seen  that  weapons  of  yield  greater  than  1 MT  ore  dictated  to 
provide  extensive  areas  of  sufficient  irradiai'on  to  rebably  ignite  window  coverings  and  room 
furnishings.  Figure  30  ^ows  peak  overpressure  versus  distance  from  ground  zero  as  a  function 
of  we^on  yield  and  detonation  altitude. 

Rgure  3 1  shows  primary,  secondary,  extinguished,  and  total  ignitions  in  tract  76  for  the 
weapems  mown  in  Figures  29  and  30.  Tlie  center  of  tract  76  was  about  5  KM  from  ground  zero. 
Thus,  for  the  1  MT  case,  an  irrudiance  of  about  4  cal/cm^,  and  greater  than  1  psi  overpre ss\^  were 
experienced.  This  accounts  for  the  moderate  number  of  primary  ignitions  ana  blast  extinctions 
sliown  In  Figure  30.  In  the  0.01  MT  case,  an  irradlancc  of  about  0.04  caiycm*’*2.  and  about  0.4 
psi  overpressure  were  experienced.  This  is  consistent  with  the  small  numbers  of  primary  ignitions 
and  blast  cx'inciions  sliown  in  Figure  30. 

It  is  sontewhat  sunnising  that  so  many  more  secondary  ignidons  than  primary  ignitions 
occur  tu  these  relauvciy  low  levels  of  overpressure.  However,  it  is  important  to  keep  in  mind  that 
ihwc  results  are  provided  for  demonstration  of  tlie  model  stiuaure,  and  can  probably  be  in  error  by 
as  much  as  an  order  of  magnitude.  A  primary  interest  in  the  debugging  process  was  to  assure  that 
the  results  exhibited  the  cottea  trends  as  funaions  of  the  imponant  variables. 

A  1 .0  MT  weapon  was  detonated  at  cooidiiiates  XB  =  19  km.  YB  =  19  km  wliii  a  burst 
height  of  HB  e  3.7  km.  Thus,  the  data  shown  in  Figures  29  and  30  apply.  The  btiundarii  j  of  the 
annulus  were  spccillcd  as  KMNUIT'  ^  20  km  and  RMXEFF  s  54  krn.  The  burst  location,  and  the 
collection  ol  all  tracts  whose  centers  are  located  in  the  unnulus  are  shown  in  Figure  32.  It  can  be 
seen  that  none  of  die  industrial/commerclol  tracts  aie  located  in  the  annulus.  This  was  done  for  two 
reasofis.  First,  it  is  considered  a  realistic  sccnaiio  for  use  of  fire  ar  a  collateral  damage  mechanism. 
Tlie  weaiwn  would  innsi  likely  be  targeted  so  that  the  desired  level  of  damage  10  the  conunerciaV 
industrial  area  would  be  caused  by  blast  effects.  It  Is  unUkclv  that  fire  would  be  depended  upon  to 
provide  sure  damage  to  those  toigeui.  Second,  the  commerclai/industrial  oacts,  as  configured, 
retjuiied  too  long  to  execute.  Tract  75,  wheii  evaluated  Individually,  riH^uired  about  50  rninutes  to 
execute  widi  duee  occupancy  classes.  Tills  long  execution  lime  can  be  avoided  by  reconfiguring 
iimi  75  iiiio  iluec  sinrUler  trucis.  each  with  one  occupancy  class.  Tluec  such  tracts  .should  execute 
in  about  1.5  minutes  (1  0.  about  0.3  muiuie  each). 

Table  i  sliows  the  list  of  tracts  Included  ht  the  annulas  of  Figure  32  (see  also  Figure  2).  It 
can  lx;  seen  diat  IV  ol  tire  44  iiucis  in  die  annulus  have  some  degr^  of  residential  occupancy,  Of 
die  19  uacis  coiuauilng  rasaleiiual  occupancy,  6  are  nredUiin  density  ( 10  to  30  jierceii;)  to  high 
density  tiU  to  lUO  ppn.ent)  otcupimcy,  while  1.3  arc  low  density  (lew  than  U)  percent)  occupancy, 


71 


y  »  yield  (MT) 

KB  »  height  of  burst  (KM) 


■y  «  1.0  MT,  HB  “  3.7  134 


y  »  0.1  MT.  HB  =  1.6  KM 


y  o  0.01  Ml’,  HB  -  0.8  KM 


Range  £ran  Ground  Zero  -  KM 


J-igure  29.  Tliemwl  Irrudiuucc  ai  dlwnnccs  from  gfoutwl 


n 


Y  s  yield  (MT) 

HB  =  height  of  burst  (KM) 


I'igure  30.  Peak  uveipresiturc  ui  ranges  from  ground  zero. 


Number  of  Ignitions 


NTOT  =  total  number  of  ignitions 
NSEC  =  number  of  secondary  ignitions 
NPRIM  B  number  of  primary  ignitions 

NEX  B  number  of  extinguished  ignitions 


Figure  1 ,  1‘ritnury,  swcoiiiiary,  and  Uilui  Ignllluii  lor  three  lieldM, 


n  BARRiH  (NO  FUEL)  ^  COME2CIAL 

□  VACANT  (NATURAL  FUEL)  ^  HIGH  RISE  RESIDENTIAL 

^  RESIDENTIAL  JU  INDUSTRIAL 


Eisu  -  Kllameuoi'i 


I'lnuru  .^2,  (Juoinelry  of  uiinulaH  uiiil  liiuluilotl  iructa, 


Tables,  Example  of  tracts  In  annulus, 


TRACT! OCCUPANCY  1 

TRACT 

OCCUPANCY 

ID 

CLASS 

ID 

CLASS 

5 

B 

60 

R 

7 

B 

87 

R 

8 

R 

88 

R 

9 

R 

09 

B 

10 

R 

90 

B 

11 

B 

91 

R 

18 

B 

105 

B 

19 

R 

107 

B 

20 

R 

109 

B 

21 

B 

110 

R 

24 

R 

113 

B 

26 

B 

114 

B 

27 

B 

115 

R 

28 

R 

120 

B 

38 

B 

121 

B 

39 

B 

128 

R 

40 

B 

131 

R 

42 

B 

132 

B 

43 

R 

133 

B 

57 

B 

134 

B 

58 

8 

59 

R 

B  -  NO  OCCUPANCY  (BARREN) 
R  =>  RESIDENTIAL  OCCUPANCY 

Y  -  1  wr,  XB  «  19  KM, 

YB  ■  19  KM,  HB  -  3.7  KM 

RNMETF  "  20  KM 
HMXEEF  «  54  KM 


SU^MARY|  44  TOACTS  IN  ANNULUS 

25  'l-RACTS  IN  ANNUUJS,  NO  OCCUPMCY 
19  TRACTS  IN  ANNUUJS,  RES.XDENTIAL  OCCUPANCY 


6  RESICENTIAL  TRACI'S  IN  ANNUUJS,  MEDIUM  TO  llicai  DEUSITI 
13  RESICENTIAL  HUOTi  IN  ANNULUS,  LOW  DENSTTY 


SECTION  6 

TRACT  BURNING  AND  PROPAGATION  MODELS 


,  ^  discussed  In  Sections  2  through  5,  the  model  of  an  urban  area  is  described  by  fuel 

cnoractetistlcs  and  geographical  location  relative  to  the  burst.  A  nuclear  burst  ignites  this  fuel  with 
Its  therrn^  pulse,  which  is  followed  by  a  blast  wave.  This  extinguishes  some  of  the  ignitions  and 
may  rubblizc  pim  of  the  fuel  bed  to  some  radius  from  Ground  Zero.  The  blast  wave  is  also 
res^rustble  for  initiating  secondary  ignitions.  ITils  is  a  statement  of  the  initial  condltloas  that  lead 
to  the  dynamic  aspect  of  tract  burning  and  the  propagation  of  fire  through  the  modeled  urban  area. 

6.1  CHEMICAL  RKACnON  MODEL. 


Tte  model  of  tract  burning  we  call  a  chemical  reaction  model,  because  its  form  is  idenUcal 
to  a  simple,  coupled  reaction  involving  three  species,  or  substances.  It  is  as  though  substance  A 
changes  into  substance  B  with  some  rate,  and  substance  B  then  changes  into  substance  C  with 
some  oUicr  rate,  'fhc  rates  at  which  these  reactions  go  can  be  dependent  on  many  things,  but  one 
ming  they  ^  certainly  dependent  on  is  the  amount  of  the  source  substances.  For  example,  if  the 
wnount  of  A  goes  to  zero,  then  B  can  no  longer  increase,  or  if  B  goes  to  zero  then  the  channel 
between  A  and  C  is  broken.  If  the  rate  at  which  B  changes  to  C  is  much  faster  than  the  rate  from  A 
to  B.  men  wc  could  simplify  by  ignorine  the  B  state  and  using  the  A  to  B  rate  as  the  effective  rate 
from  A  to  C.  Ihe  numerical  differential  equations  that  describe  this  system  arc  integrated  forward 
m  time,  accounting  for  changing  conditions  that  affect  the  rates. 

In  the  Uact  bvuiiing  model  the  three  quantities  are.  for  tract  Index  1, 


the  number  of  urtlgnitcd  buildings 


the  number  of  burning  buildings 


the  number  of  burnt  out  buildings. 


n‘i 

Ni 


llieir  uniform  ateal  deasitles  are 
n'l  s  Nj  /  Ai 

ni  s  Njj  /  Ai 
n3  s  /  Ai 


(33) 


where  Ls  urcu  of  tract  1  in  square  meters.  Note  iliui  the  cunscrvmiun  of  building, i  gives 


N'  ^  n‘i  4.  Ni  4-  Ni 


(34) 
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where  Is  the  total  number  of  buildings  in  the  tract.  Now  the  initial  conditions,  with  primary  and 

secondary  ignitions,  arc  that  ^3  "  ®  and  that  ^1  “  ^  *  ^2c ,  where  represents  the  initial 
number  of  ignited  buildings. 


To  write  down  the  governing  equation  for  ^1,  which  involves  the  propagation  of  fire  to 
unbumed  structures,  we  equate  the  rate  at  which  unbumed  structures  are  ignited  to  the  product  of 
their  effective,  average  cross  section;  the  number  of  unbumed  buildinp;  and  the  rate  at  which 
brands  that  would  cause  ignitions  will  land  per  unit  area.  Symbolically,  this  becomes 


dN‘i 

dt 


-ai  n‘i 


All  1  y  Aid  1 

Ai  t^Ai  k*i  't'^Aic 


or 


ciN‘, 

dt 


-oi  n‘j 


^1  +  y 
Ait^  h 


Aid  1 

Ale 


(35) 


(36) 


where 

is  the  effective,  average  cross  section  of  each  building  in  tract  i 
All  is  the  arca  of  tract  i  that  is  die  source  of  brands  that  land  in  track  i  suming  fires 
Aw  is  the  area  of  dact  k  that  is  the  source  of  brands  Uiat  land  in  track  i  suning  fires 
Is  the  total  brand  production  period  for  tract  1, 


Tljus  All  /  Ai  represents  the  fraction  of  duct  i  urea  that  contributes  to  new  fires  in  duct  i.  Akl  f  Afc 

rcprescnis  the  fraction  of  dact  k  uteu  that  contributes  to  new  fires  in  tract  i,  ^  ^  Ai) 

1 

brand  production  rate  per  unit  area  in  duct  i  over  all  occupancy  classes,  and  is  the  probability 

Uiut  u  brand  landing  in  will  start  u  new  fire.  Only  ducts  Uiat  have  a  common  boundary  are 
considered  for  conUigion,  which  is  represented  in  the  sum  to  three.  Obviously  tracts  on  the 
boundary  of  the  urban  area  modeled  may  have  two  or  only  one  neighboring  tract,  so  the  sum  to 
tiuee  is  symbolic  of  the  typical  cu.se,  and  is  not  to  be  iniei'prcied  as  a  sum  over  the  first  three  tracts. 

Tlic  other  two  govcmlng  equations  arc 


dNi  dN 


dl 


dt 


w,' 

I  d7 


evaluated  at 


t  -  t;, 


and 


(37) 


dt  \  dt  / 


(3K) 


7H 


evaluated  at 


T 

where  ^  o  is  the  average  burnout  time  for  buildings  in  tract  i,  which  is  the  weighted  average  over 
occupancy  classes.  In  this  version  of  the  Urban  Fire  Simulation  (UFS-2)  the  burnout  time  for  a 

b'iilding  in  occupancy  class  j  is  a  constant  ^  o.  Tlius  the  rate  at  which  new  class  j  fires  started  in 

U'acl  i  0  seconds  ago  is  also  the  rate  at  which  hvaming  buildings  are  now  reaching  the  burnt  out 
state.  Obviously  the  model  could  be  enhanced  at  this  point  for  added  realism,  but  as  it  stands  the 
UFS  -2  is  adequate  for  demonstrating  the  process  from  ignition  to  final  burnout.  One  such 
extension  would  be  to  account  for  multiple  ignitions  per  structure,  which  would  lead  to  a  quicker 

burnout  and  a  more  rapid  increase  in  ^o,  the  heat  production  rate. 

The  branding  rate  per  building,  or  its  reciprocal  the  branding  period  in  seconds,  is  taken  as 
a  function  of  the  average  heat  production  rate  in  the  burning  building  given  in  watts 


Qo 

iBo  (39) 

This  branding  rate  is  100  times  the  rate  used  in  die  previous  Urban  Fire  Simulation  (Ref.  1),  This 
is  consistent,  because  the  average  branding  range  a.ssumcd  in  the  previous  work  was  about  1  km. 
because  of  the  grid  cell  size.  In  the  model  described  here  a  distance  of  1(X)  m  is  much  more 
consistent  with  Llie  concepts  modeled.  'I'o  extend  tliis  branding  riue  to  the  total  branding  rate  for  a 
tract.  SCI 


I  ub  classes 


lOM  X  n!/Qo  =  IO'^/Qt 


where  Q 1’  is  the  total  heat  production  rate  for  tract  i. 

Hie  effective,  average  building  cross  section  is  computed  as  follows: 


all  clusses 


I  N'i  / 1,1 


svlicrc  fi)  is  the  CratUonaJ  area  of  tract  i  covered  by  buildings  ol  clu.ss  j,  I  lowever.  In  the 
compulations  in  UFS-2  ibe  lenii  ^1  1  /  is  developed  by  evaluuling 


.,1  ail  classes  all  classes  .,i| 

A|  I  J  N'J 

As  u  simple  lesi  ol  the  iuuiidlng  rule  eonsluiil  an  urban  icsidenllul  liucl  ( I  ‘KKMI  sltuelurch  In  20  kiii' )  with 
It  «  10(111  nv’,  was  iiilllully  .M)  peieeni  luiulondy  Ignlieil,  In  24  mhiules  ‘hi  |iereenl  ol  llie  MiucUites 
III  the  liuu'l  wcie  nil  llie. 


For  each  of  the  initial  fire  starts,  priniary  and  secondary  ignitions,  a  burnout  event  is 
•T*j 

scheduled  at  o.  Similarly,  ignitions  that  are  the  result  of  branding  at  time  t  will  have  burnout 

t  4- 

events  scheduled  at  o  for  each  building.  Another  obvious  extension  to  UFM  would  be  to 


draw  0  from  a  distribution  rather  than  treat  it  as  a  constant,  or  to  update  the  burnout  fraction  of 
each  burning  building,  and  switch  it  to  category  3  when  this  fraction  exceeded  .99.  This  last 


method  would  allow,  for  example,  wind  to  affect 


rpij 

as  well  as  0. 


Wind  in  UFM  is  the  vector  sum  of  ambient  wind,  which  the  user  defines,  and  winds 
induced  by  the  fire.  Conceptually  a  tract  would  have  an  image  associated  with  it  (for  a  short  while) 
in  the  moving  air  mass  above  it.  Refer  to  Figure  33,  where  the  displacement  is  greatly 

exaggerated.  The  overlapped  areas  are  proportional  to  the  rate  of  self-branding,  tract  i  into  tract 

i,  and  contagion,  tract  k  into  tract  i  and  ■^'j  tract  i  into  tract  j. 

When  updating  tract  i,  the  branding  rates  must  be  computed.  Because  the  wind  velocity 
differs  from  tract  center  to  tract  center,  using  the  wind  at  the  center  of  tract  i  to  estimate  the 
branding  that  occurs  at  tract  i  boundaries  would  introduce  errors.  If  only  the  wind  at  tract  i  is 
considered,  then  the  contagion  branding  rates  will  not  agree  with  those  computed  when 
neighboring  tracts  arc  updated.  For  example,  assume  that  two  tracts  like  i  and  j  in  Figure  33  lie 
west  and  cas»  respectively,  and  that  the  wind  in  tract  i  is  blowing  to  the  southeast  while  the  wind  in 
tract  j  is  southwest  In  sucl>  a  case  an  update  of  tract  i  would  show  brands  landing  in  tract  j  and 
vise  versa.  A  better  solution  is  to  make  a  lincai  estimate  of  the  wind  at  the  ij  boundary  and  use  this 
wind  to  piedict  the  branding  along  that  boundary.  This  is  in  fact  what  is  done  in  UFS-2.  The 
overlapped  areas  arc  computed  a.s  tlie  product  of  the  length  of  a  boundary  and  the  component  of 
wind  perpendicular  to  the  boundary  for  contagion.  Self-branding  then  ticcomes 

3 

Aii  =  Aj  '  Aki  +  Aik 

k*i  (43) 


wiicre  cither  Aki  or  Aik 


will  be  zero  for  each  k. 


I'lguic  33,  Ulspluccineiu  of  iiir  muss  Ixrcuusc  of  wlitds, 


HU 


Contagion  is  assumed  to  occur  just  along  the  boundary  where  branding  exists;  it  is  not 
averaged  over  the  tract  area.  This  would  result  in  a  serious  numerical  diffusion  error  in  simulated 
fire  propagation.  The  spread  of  fire  across  a  tract  that  had  been  previously  unignited  would  be 
effectively  instantaneous.  Because  this  would  seriously  violate  the  intent  of  UFS-2.  an  average 
propagation  time  is  kept  and  updated  for  the  progress  of  the  fire  front  across  a  tract  ignited  by 

contagion.  This  propagation  lime  delay  ^  ^  P  for  tract  j  ignited  from  tract  i  is 

ATp«  =  ^  va; 

(44) 

where 

Lt  -  Vw  ffii  Ls  the  effective  branding  range  (m) 

is  the  wind  velocity  perpendicular  to  the  ij  boundary  (m/s) 

I  .  i  i  \l/3 

tai  -0.1  (Qx  /  N2I  ^  average  brand  flight  time  (s)  (45) 

and 

A.ij  =  Li  Sjj^  branding  range  limes  the  ij  boundary  length  (m2), 


When  the  fire  front  has  crossed  such  a  tract  then  it  can  centribute  brarids  uO  its  neighbors. 

An  Allowed  Propagation  from  tract  j  event  is  scheduled  at  where  is  die  current  update 

lime.  Subsequent  tract  updates  will  include  more  branding  from  tract  i,  and  perhaps  other  tracts  as 
well,  and  self-branding  of  tract  j.  Because  the  branding  rates  and  flight  times  are  functions  of  the 
rate  of  energy  release,  a  fire  that  Is  increasing  rapidly  in  intensity  can  cause  the  fire  front  to 
accelerate  across  a  tract  To  account  for  this  a  tract  that  is  prevented  from  propagating  the  fire  to  its 
neighbors  will  have  Allowed  Propagation  events  stored  at  each  update.  The  first  of  these  that  is 
encountered  during  the  simulation  will  allow  the  mact  to  conuibuie  to  fire  spread;  all  others  will  be 
ignored. 


To  carry  out  a  u-act  update  the  generating  differential  equations  must  be  UTinsformcd  into 

their  numerical  equivalents.  The  equation  for  ^3  simply  becomes  a  tallying  of  burnout  events. 

The  equation  for  ^2  is  the  accumulated  result  of  the  new  fire  sums  since  the  last  update  reduced  by 

the  number  of  burnouts.  And  the  equation  for  the  change  in  is  the  product  of  the  average 
ignition  rate  and  the  update  interval.  Symbolically 


AN  =  -A 


tiN*. 


2  cll 


dN; 

“dT 


U) 


(l  - 1  p 

\  II  ml  ^ 


(46) 


where  the  update  time  has  udvimeed  from  to  P  is  u  uniform  random  number  on  (0,1),  and 

^^1,  which  i,s  defined  as  a  positive  i)un''v is  cun.sirained  by  ^1,  q'jn;  random  number 

(3  causes  an  appropriate  rounding  to  whole  numbers  of  buildings.  I'or  example,  3,2  fire  .starts  in 


»1 


the  update  intei’val  will  round  to  4  furc  starts  20  percent  of  the  tinic  and  to  3  fire  stais  80  percent  of 
the  lime.  Tlius,  at  update  time  is  reduced  by  ^2  is  increased  by  and  reduced  by 
the  number  of  bum  out  events  in  the  update  interval,  and  just  accumulates  bum  out  events, 


The  next  update  time  for  tract  i  is  estimated  by  linear  extrapolation  as  the  time  by  which  a 

10  percent  change  could  occur  in  ^1,  ^2,  or  ^3  with  a  minimum  and  maximum  allowed  lime 
step.  The  actual  fonnulation  that  appears  in  LIF.S'2  is  as  follows; 


dNl 


=  0 


implies 

Ni  =  0 


N‘ 

at;  =0.1  -  * 


=  600  dN;  /  dt 


=  600  at;"  =  min  [600.  6  N‘  /  N; 


AT‘  =  min  {l800,  max  (sO,  AT;  ,  AT;  ))  (47^ 

Note  that  tracts  are  not  updated  in  lockstep  at  some  sp^ified  tnterval.  but  are  updated  individut^y, 
asynchronously.  The  next  update  event  for  each  burning  tract  i.s  placed  on  the  event  list  at  the  lime 
of  the  current  update.  ITiis  asynchronous  updating  is  a  fcatuiu  (if  the  UFS-2  design.  It  provides  a 
combination  of  good  accuracy  and  computational  efficiency.  Tne  only  exception  to  this  is  that  all 
tracts  are  updated  prior  to  an  output  event  reque.stcd  by  the  user. 


The  energy,  power,  and  bum  out  time  per  building  in  each  occupancy  class  is  given  in 
Table  4.  The  "Vacant"  class,  which  represents  naiurul  fuel,  has  been  left  blank  in  the  current 
model.  It  could  be  filled  in.  for  example,  with  tm  appropriate  energy  per  tree,  power  per  tree,  and 
bum  out  time  for  a  uee  and  the  number  of  trees  in  a  vacant  tract  v.'ould  replace  the  number  of 
buildings  on  input. 


Blast  from  the  nuclear  weapon  can  rubblize  part  ot  ilie  fuel  in  the  urban  area.  If  a  tract  is 

qd  fV  • 

i'ubbiized,  the  effective  number  of  buildings  in  each  class  is  halved,  o  is  halved,  and  is 

reduced  to  one  third  of  its  nonnol  value,  llius,  ^  is  one  sixth  of  the  energy  [Tcr  building,  which 
allows  only  one  iwelfUi  of  the  louil  energy  to  be  released  when  the  uact  is  burnt  out. 
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Table  4.  Approximate  constants  for  rectangular  bum  curves  by  occupancy  class. 


Occupancy 

Class 

Qo  (joules) 

()o  (watts) 

To  (s) 

Barren 

0 

0 

0 

Vacant 

0 

0 

0 

Residential 

. . . . 

2(10)^° 

4(10)^ 

5(10)3 

Commercial 

8(10)‘° 

8(10'/ 

KIO)" 

- ntginuis — 

Residential 

1(10)‘‘ 

1(10)’ 

KIO)'* 

Indusuial 

6(10)^^ 

6(10)’ 

1(10)“^ 

6.2  FIRE  INDUCED  WIND  MODEL. 

Surface  winds  resulting  from  a  fire  column  are  an  inflow  response  to  the  enuairunent  and 
lofting  of  surface  aii.  In  the  UFS-2  wind  model  each  burning  tract  is  considered  to  be  a  subfire 
with  its  own  fire  column.  Nieghboring  subfire  columns  coalesce  at  some  ^titude  above  the 
ground,  and  this  larger,  combined  column  spreads  in  a  conical  fashion  as  it  rises.  The  center  of 

J 

each  triangular  tract  is  regarded  as  the  subfire  center,  with  an  effective  radius  f  given  by 
Tf  -  (48) 


Each  subfire  is  made  up  of  ihe  fire  columns  from  each  of  the  burning  buildings  in  the  unict,  fm* 
columns.  The  coalescence  height  of  fire  columns  u: 

4  ==  4.6263 

\)tN2/  (49) 

where  (here  is  more  than  one  column,  and  the  origin  displacement  below  the  surface  is 


4  =  4.6263  4  -  4 

in  meters,  Iliese  disuuices  are  constrained  by 

0  <  zj,  <  4.6263  r| 


(30) 


(51) 
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and 


0  <  Zc  <  4.6263  r[  (-52) 

Geometrically,  the  coalesced  fire  columns  in  a  subfire  are  depicted  in  Figure  34,  In  the 

construction  of  Figure  34  about  ten  tire  columns  were  assumed,  which  makes  ^  about  twice 
In  addition,  to  make  uhe  figure  more  legible  the  horizontal  scale  has  been  exaggerated  by  about  two 
times. 


Figure  34.  Geometry  of  coalescing  fire  columns. 


The  radial  indrali  component  of  .surface  wind  is  defined  us  a  vector  pointing  to  llic  .subfne 

(tract)  center.  Let  these  locution  vectors  be  given  by  and  ilie  test  point  where  indriU't  is  to  lx; 

evaluated  given  by  ^r.  Tlic  radial  component  is  directed  olon^  the  vector  "  ^^i'.  .\l'icr  a  wind 
velocity  is  assigned  to  each  of  these  vectors  lone  for  each  subfire)  they  are  .summed  to  obtain  die 
combined  indraft  wind.  Let 


r  = 


•  Rp 


be  die  distance  from  the  center  of  tract  i  to  the  test  point. 


(33) 


K4 


If  ^  the  test  point  lies  outside  the  tract  i  subfire,  then  the  normalized  velocity  of  the 
indraft  wind  is  given  by 


vkr) 


Va 


3  ak 

/  •  i  \ 

2i 

5  Ti 

1/3 


i  Tf  i2J3 

Ni4zc  + 


^  2i  [f(4,(o‘)  -  f((o>')] 
T1  \  r  / 


(54) 


where  the  first  term  is  the  air  required  by  the  burning  rate,  the  second  term  is  the  air  entrained  by 
the  lower  column,  and  U.e  third  term  is  the  air  entrained  by  the  upper  column.  The  new  symbols  in 
Equation  54  are  as  follows: 


V3  =  ^  ^ 

Q  ^  is  the  volume  rate  of  air  required  (m3/s)  (55) 

where 

A  is  the  mass  of  air  needed  to  bum  a  mass  of  fuel  (kg  air  /  kg  fuel) 

Q  is  the  heat  produced  by  a  mass  of  fuel  (J/kg) 

d  is  the  air  density  (kg/m3) 

ot  —  0. 1 ,  k  =  3.897(10)  and  T)  =  5.5 1  constants 
and 


ro)T 


ii2/3 


0 


(a  ■*-  MqI 

(1  + 


dco 


(56) 


is  the  truncated  column  entrainment  factor,  a  tabulated  function,  where  the  scaled  reciprocal 
distances  are 


(Dj  —  z-p  /  r 
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(57) 


and 


I  .  j  j  3/4 

Zj  -  0.234  Qx  ATq  jg  column  top  from  tract  i  (m)  (58) 

^  B  is  the  length  of  time  tract  i  has  been  burning. 

The  UFM-2  model  enforces  the  following  constraints  on  the  truncated  column  enirainment 
parameters: 

0  ^  oDq  5  4.6263  ^59) 

0  <  CDc  <  COt  <  50 

A  plot  of  is  shown  in  Figure  35. 


•  1  Z  3  4  S  » 

Scaled  rec^rocal  dtetance  to  coKinn  origin  -w. 

Figure  35.  Truncated  column  enu  ainment  function. 

If  ^  then  the  test  point  is  inside  the  tract  i  subfire  and  the  wind  velocity  from  this 
subfire  is  scaled  linearly  from  the  tract  center.  The  wind  •.■elocity  is  given  by  Equation  54 
i 

evaluated  at  *'<'  times  the  distance  to  the  tract  center; 
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(60) 


vi(r)  =  vkrj) 


The  radial  indraft  and  the  ambient  wind  is  then  summed  over  all  velocity  components  to  obtain  the 
wind  at  test  point  p 


ail  burning 
tracts 

v;(r‘) 

r  (61) 

which  can  never  have  a  zero  distance  for  ^  when  the  lest  point  is  on  a  tract  boundary.  The 
contagion  area  can  then  be  expressed  as 

Aii  =  tfiil  Vpx  .SijI  (62) 

where  is  a  tract  edge  vector. 


6.3  UFM-2  OUTPUT  EXAMPLES. 

As  a  final  test  of  the  UFM-2  software,  three  test  problems  were  executed.  These  tests  used 
^e  urban  map  and  tract  layout  shown  as  an  example  in  Sections  2  through  5.  The  three  tests 
involved  three  different  ambient  wind  conditions:  no  wind,  1  m/s  wind,  and  10  ra/s  wind  from  the 
southeasL  A  1 .0  mT  burst  is  assumed  to  occur  over  the  boundary  of  tracts  52  and  77 
( 19km,  19km)  at  an  altitude  of  3700  meters.  All  three  tests  completed  successfully  with  expected 
results. 


Vp-v.-b 

i 


Although  almost  all  of  the  occupied  tracts  are  on  fire  by  2000  seconds,  the  major 
conuibutors  to  energy  release  rate  are  uhe  residential  tracts  33,  34,  and  50,  and  the  high  rise 
residential  and  indusuial  tracts  67,  68,  69, 100,  and  101.  There  is  little  difference  between  the 
zero  ambient  wind  case  and  the  1  m/s  wind  case.  Tract  bum  out  occurs  in  the  same  order  an  with 
insignificant  differences  in  burnout  times,  which  results  from  the  Monte  Carlo  simulation  design. 
The  last  five  tracts  listed  above  are  also  the  first  to  bum  out;  u^ct  68  at  about  5400  seconds  is  the 
first  The  only  apparent  difference  in  these  two  cases  is  that  tract  108  caught  fire  from  tract  109 
brands  at  6540  seconds  presumably  because  of  the  1  m/s  southeast  wind.  This  wind  component  is 
in  the  right  direction  to  influence  that  contagion,  while  radial  inflow  alone  would  not.  Tracts  108 
and  109  are  in  the  southwest  comer  of  the  urban  map.  Bear  in  mind,  however,  that  the  Monte 
Carlo  simulation  design  would  allow  this  contagion  event  to  occur  earlier,  later,  or  not  at  all  in 
other  trial  runs  using  different  pseudo  random  number  seeds.  The  major  reason  for  the  similarity 
of  results  is  that  the  largest  indraft  component  at  a  tract  center  at  2000  seconds  is  13  m/s  and  is 
10.5  m.s  at  4000  seconds.  This  indraft  clearly  doramaies  in  the  principal  fire  areas. 

In  the  10  m/s  ambient  wind  case  the  order  of  tract  burnout  is  affected  to  some  degree,  and 
the  ambient  wind  sped  up  the  burning.  The  first  tract  to  bum  out  is  138  at  5140  seconds  followed 
by  tract  68  at  5250  seconds.  Tracts  that  quit  burning  early  in  the  simulation  typically  leave  ten  or 
fewer  suuctures  unbumed  out  of  several  thousand.  Tracts  that  quit  burning  toward  the  end  of  the 
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simulation  (18000  seconds.  5  hours,  in  this  case)  typically  leave  a  few  percent  to  about  twenty 
percent  of  structures  unbumed. 

Output  from  the  1  m/s  ambient  wind  case,  that  is  plots  of  quantities  recorded  at  requested 

simulation  times  and  tract  maps  plotted  by  UFM'2  of  relative  wind,  tract  burning,  and  Q,  are 
shown  in  Figures  36  through  53.  Figure  36  is  a  plot  of  the  number  of  tracts  unignited,  burning, 
and  burned  oul  Figure  37  is  a  plot  of  the  maximum  wind  speed  at  tract  centers.  Figures  38 
through  42  arc  tract  maps  of  relative  wind  vectors.  Figures  43  through  47  are  uact  maps  of  tract 

burning  and  burned  out  status.  Figures  48  through  53  are  uact  maps  of  Q.  The  values  of  Q 

associated  with  each  tract  arc  coded  on  a  logarithmic  scale,  base  10  with  an  offset  of  .  Thus,  1 
means  l.e5,  2  means  i.e6,  etc. 

Other  parametric  studies  with  UFS-2  involved  multiple  detonations  and  patterns  of  yield 
and  altitude,  wind,  and  fuel  loading  that  might  dis''lose  modeling  errors.  One  such  pattern  placed 
the  bursts  over  the  northeast  comer  of  the  urban  map  and  used  single  detonations  of  0.(K11  mT  at 
325m.  0.01  mT  at  750m,  0.1  mT  at  1500ra  and  l.o  raT  at  3500m.  Another  pattern  was  a  2x2 
matrix  of  low  and  high  ambient  wind  and  low  and  high  fuel  loading.  All  of  these  tests  were 
performed  satisfactorily. 
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Figure  37.  Maximum  wind  speed  at  tract  centers. 
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Figure  38.  Tract  map  of  relative  wind  vectors  at  2000  s. 
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Figure  39.  Tract  map  of  relative  wind  vectors  at  4000  s. 
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Figure  40.  Tract  map  of  relative  wind  vectors  at  6000  s. 
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Rgurc  41.  Tract  map  of  relative  wind  vectors  at  8000  s. 
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Figure  42.  Tract  map  of  relative  wind  vectors  at  10000  s. 
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Figure  43.  Traci  map  of  unignited,  burning,  and  burned  out  tracts  at  2000  s. 


95 


T  AXIS  (KlLOMETERSl 

.57  H.I«  16.71  22.29  27.06  3.3.113  39.  DO 


HISSIQN  RE5ERRCH  CQRPORPT  J 

URBAN  FIRE  MODEL  VERSION  2 


PLOT  SlUUUATlON  TIME  -  4000.00  SECONDS 

nw  smt  MOft  B««fWiNC  tract.  x«mr.NCD  cut,  oncRKisc  uhiciiitcd 
HinBcn  or  uNiCHirco.  burmihc  swi  Buwto-eur  tracts  «  86  rs  c 


"^-OD  5.43  l'0.B6  r6.29  21.“:  27.14  32.57 

X  RXI5  IKILOMETEBSj 


Figure  44.  Tract  map  of  unigniied.  burning,  and  burned  out  tracts  at  4000 
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Figure  45.  Tract  map  of  unigniied.  burning,  and  burned  out  uracts  at  6000  s. 
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Figure  46.  Traci  map  of  unignited,  burning,  and  burned  out  tracts  at  12000  s. 
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Figure  47.  Tract  map  of  unigniied,  burning,  and  burned  out  tracts  at  24000  s. 
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Figure  49.  Tract  map  of  Q  at  4000  s. 
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Figure  50.  Tract  map  of  Q  at  6000  s. 
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